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INTRODUCTION acceptance of his conclusions should be withheld. 
Rie Sip aap ; Since a clay effect is possible, Odeh should have 
MS The general application of Darcy’s law to natural been allowed space to present a more detailed 
a rocks has already been challenged in the literature. ! account of the rocks which he used and the pre- 
mart The evidence shows that the permeability as calcu cautions which he took to avoid confounding the 
lated from the Darcy equation can be a function of effects of clay and fluid circulation. In addition, 
the pressure drop and the salt concentration of the it follows from Odeh’s discussion that the oil 
water phase. Most explanations for aberrant behavior relative permeability should increase as the vis- 
involve clays and their properties and have wen cosity ratio increases, Odeh presents data corrob- 
" qualitatively satisfactory.Recently, however, Odeh orating this deduction, It also follows that the oil 
en- revived the theoretical views of Yuster.? Since relative permeability should decrease as the vis- 
At- Yuster’s concept implies a fundamental error in cosity ratio of oil to water decreases. However, our 
is using Darcy’s relationship for two-phase flow, and results show that a high relative permeability, when 
at not merely that conditions may limit its use, Odeh’s once attained by using a viscous oil, may be main- 
“a experimental support arouses considerable interest. tained when that oil is replaced by an oil of much 
‘“ However, Odeh’s work as presented is thought to lower viscosity, 
al to be inadequate. He has omitted important informa- 
pir tion about his materials and procedures; therefore, EXPERIMENTAL PROCEDURE 
a- 
ni- Original manuscript received in Society of Petroleum Engineers The results are shown in Fig. 1 and were obtained 
d- office Sept. 27, 1960. Revised manuscript received Nov. 30, as follows. Sandstone cores, 3-in long x l-in. in diam- 
1960. ; ; . , 
ht IReferences given at end of paper. eter and which contained clays, were saturated with 
m 
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FIG. 1 —-— EFFECTIVE PERMEABILITY RELATED TO THE SEQUENCE OF OPERATIONS. 
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brine (50g. NaCl/l. solution),and the permeabilities 
were measured in a slightly modified Penn State 
apparatus. Each core was flooded with the lowest- 
viscosity oil (i.e., Soltrol having a viscosity of 
1.42 cp), and the oil permeability was measured 
with the residual brine stationary. This was the 
saturation region which Odeh found produced the 
greatest effect. The low-viscosity oil was then 
displaced by higher-viscosity oils — either Wemco 
(viscosity = 13.2 cp) or Merusol (viscosity =73.3 cp) 
— and the effective permeability to oil was again 
measured. The low-viscosity Soltrol was again 
introduced into the core so as to displace the Merusol 
miscibly, and the effective permeability was meas- 
ured, The fluids followed one another in the order 
shown in Fig. 1. 


CONCLUSIONS 


We conclude that oil viscosity can influence the 


effective permeability of some natural rocks to oil, 
However, when once attained, such an increase 
is not necessarily lost on replacing the high-vis. 
cosity oil with one of low viscosity ~ an observation 
strongly at variance with Odeh’s analysis but 
qualitatively explicable in terms of the movement 
of colloid particles at oil-water interfaces. 
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DISCUSSION 


A. S. ODEH 
MEMBER AIME 


It is gratifying to see that Downie and Crane, on 
the basis of experimental results, concluded that 
‘‘oil viscosity can influence the effective permea- 
bility of some natural rocks to oil’. The authors 
also suggest that once the influence is attained it 
is not lost, which is at strong variance with Odeh’s 
analysis. 1 

In examining the authors’ results, Odeh does not 
see such a strong variance; in fact, he detects 
support. Admittedly, in Cores 1, 2 and 3 the effec- 


AUTHORS’ REPLY 


Odek1 omitted to describe what would seem an 
obvious and necessary extension of his investiga- 
tion of the relationship between relative permeability 
and viscosity ratio. If Odeh’s analysis is correct, 
the relative permeability to oil must increase when 
the viscosity ratio increases and must decrease 
again by the same amount when the viscosity ratio 
is returned to its original value. 


However, our results indicate that this is not 
necessarily so. We presented a group of results 
(Fig. 1) which provide examples of slight decreases 


in effective permeability (Cores 1 and 2), no change 
in effective permeability (Core 3) and an increase 
in effective permeability (Core 4) when a more- 
viscous oil is replaced by a less-viscous oil. 





lodeh, A. S.: ‘‘Effect of Viscosity Ratio on Relative Permea- 
bility’’, Trans., AIME (1959) Vol. 216, 346. 
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tive permeability to Soltrol was higher in the second 
cycle than in the first, but it was still lower than 
the effective permeability to Merusol, the more 
viscous oil. This is in line with Odeh’s analysis.! 


The writet has no explanation for Core 4. 


As far as the effect of clay on Odeh’s work, Table 
11 indicates that the corrected air permeabilities 
checked with the liquid permeabilities, thus imply- 
ing a negligible clay effect if any. 


TO A. S. ODEH 


Our contention, which is in danger of misrepre- 
sentation by Odeh, is that the increase in permea- 
bility is not necessarily lost on replacing the high- 
viscosity oil with one of low viscosity. This is at 
variance with Odeh’s analysis. Support by selection, 
as employed by Odeh in his discussion of our paper, 
seems an approach of contestable validity. 

The effect of the presence of clay particles on 
the permeability of a rock is still incompletely 
understood. The identity of air and liquid permea- 
bilities is open to various interpretations until the 
liquid is more fully described. For example, if the 


liquid used were distilled water the implication 
would be that clays were absent or inert and, if the 
liquid were a concentrated CaCl, solution, that 
clays might or might not be present. However, at 
this stage in our knowledge of clays, there is no 
implication that clays (if present) would not affect 
the relative permeability to oil. lanes 
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An Imbibition Model—Its Application to Flow Behavior 
and the Prediction of Oil Recovery 


J. NAAR 
JUNIOR MEMBER AIME 


J, H, HENDERSON 
MEMBER AIME 


INTRODUCTION 


The displacement of a wetting fluid from a porous 
medium by a non-wetting fluid (drainage) is now 
reasonably well understood. A complete explana- 
tion has yet to be found for the analogous case of 
a wetting fluid being spontaneously imbibed and 
the non-wetting phase displaced (imbibition). During 
the displacement of oil or gas by water in a water- 
wet sand, the porous medium ordinarily imbibes 
water. The amount of oil recovered, the cost of 
recovery and the production history seem then to 
be controlled mainly by pore geometry. The influence 
of pore geometry is reflected in drainage and imbibi- 
tion capillary-pressure curves and relative permea- 
bility curves. Relative permeability curves for a 
particular consolidated sand show that at any 
given saturation the permeability to oil durin 
imbibition is smaller than during drainage.) 
Low imbibition permeabilities suggest that the 
non-wetting phase, oil or gas, is gradually trapped 
by the advancing water. 

This paper describes a mathematical image 
(model) of consolidated porous rock based on the 
concept of the trapping of the non-wetting phase 
during the imbibition process. The following items 
have been derived from the model. 

1. A direct relation between the relative per- 
meability characteristics during imbibition and those 
observed during drainage. 

2. A theoretical limit for the fractional amount 
of oil or gas recoverable by imbibition. 

3. An expression for the resistivity index which 
can be used in connection with the formula for 
wetting-phase relative permeability to check the 
consistency of the model. 

4. The limits of flow performance for a givea 


rock dictated by complete wetting by either oil 
or water. 


5, The factors controlling oil recovery by imbibi- 


tion in the presence of free gas. 





Original manuscript received in Society of Petroleum Engineers 
office July 13, 1960, Revised manuscript received Jan, 20, 1961. 
Paper presented at 35th Annual Fall Meeting of SPE, Oct. 2-5, 
1960, in Denver, 
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The complexity of a porous medium is such that 
drastic simplifications must be introduced to obtain 
a model amenable to mathematical treatment. Many 
parameters have been introduced by others in 
“‘progressing’’ from the parallel-capillary model 
to the randomly interconnected capillary models 
independently proposed by Wyllie and Gardner3 
and Marshall4. To these a further complication 
must be added since an imbibition model must trap 
part of the non-wetting phase during imbibition of 
the wetting phase. Like so many of the previously 
introduced complications, this fluid-block was 
introduced to make the model performance fit the 
observed imbibition flow behavior. 


THE IMBIBITION MODEL 
ASSUMPTIONS 


The following assumptions are common to many 
two-phase models previously described. 

1. The wetting characteristics of the medium are 
clearly defined. 

2. When two fluids co-exist in the medium, the 
wetting phase will fill the smaller pores and the 
non-wetting fluid the bigger ones. 

3. The pores have a circular cross section. 

4. Apparent pore radii are related to the capillary 
pressure (drainage) by 


(1) 


5. The size distribution of the pores defined 
in Assumption 1 can be derived from the drainage 
capillary-pressure curve. 

A model valid for the imbibition process must 
meet the following additional requirements. 

1. When the medium has been desaturated to its 
irreducible wetting-phase saturation S,,; , the relative 
permeability to the non-wetting phase must be the 
same whether computed for the drainage process of 
for the imbibition process. 

2. A mechanism must be provided for the wetting 
phase to by-pass pores full of non-wetting fluid. 

3. Valid results must be given when other prop- 
erties of the medium are computed. 





DESCRIPTION 


The imbjbition model is a combination of Wyllie 
and Gardner’s3 and Marshall’s4 representation of 
a porous medium, with additional provisions to 
allow the invading fluid to by-pass and block pores 
full of non-wetting phase. 

The random interconnection of pores, characteristic 
of a porous material, is created by slicing the model 
normal to the capillaries as shown in Fig. 1. Two 
parts (1 and 2) thus are created with exposed Faces 
A and B, which are rejoined to form one element 
of the model after a random rotation of Parts 1 and 
2 about the axis of the capillaries. A fluid particle 
entering the element through a capillary of Part 1 
is dispersed through pores of Part 2 at the contact 
of A and B. 

The ratio @ of void area to total area in Faces 


A and B is taken equal to the porosity of the teal 
medium. This is justified by the fact that, in an 
isotropic porous medium, space and surface porosity 
are equal. The irreducible wetting-phase saturation 
is assumed to be part of the solid matrix, and the 
porosity @ is replaced by an effective porosity ¢* 
equal to ¢ (1 — S,,;). 

The effective pore-size distribution of the model 
is defined by the drainage capillary-pressure curve 
of the real medium. To this point the model is like 
that described by Wyllie and Gardner 

The blocking of the non-wetting phase during 
imbibition can be readily visualized for a real 
medium. Access to a large pore such as that of Fig. 
1 (a) may be controlled by small pores. When all the 
adjoining small pores are filled with the wetting 
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FIG. 1—SCHEMATIC DIAGRAM OF MODEL COMPO- 
NENTS, THEIR SATURATIONS DURING DRAIN- 
AGE AND IMBIBITION. 


phase, the non-wetting fluid occupying the large 
pore cannot escape by any simple mechanism ip- 
volving viscous flow. Forms of mass transfer such 
as diffusion of gas through the wetting phase are 
excluded from consideration. The blocking mechanism 
which occurs in the real medium cannot be duplicated 
in a single element of the model such as shown in 
Fig. 1 (b). A further difficulty arises when two 
identical elements are placed in series (Fig. 1c), 
The permeability of the resulting element is not 
equal to that of the individual components, as in 
the case of a real medium. 

To overcome the latter deficiency and to satisfy 
the blocking requirement, an artificial space is 
introduced at the end of each element which will 
be called the ‘‘collecting space’’ (see Fig. 1d). 
The collecting space does not represent a physical 


teality but is used only to meet the model requires 


ments specified earlier in the paper. 

The collecting space introduces the fluids, leav- 
ing Part 2 of any element into the proper capillaries 
of Part 1 of the following element. Thus, all the 
elements are in the same condition. In addition, 
the collecting space acts as a valve to trap the 
non-wetting phase in the manner explained as 
follows. 

During the drainage process, when the capillary 
pressure becomes sufficiently large that the non- 
wetting phase invades capillaries with radius 1, 
all larger capillaries are desaturated and there is 
no by-passing of blocking of the fluids. During the 
imbibition process, the collecting space allows 
all the capillaries of Part 1 of any element to be 
saturated by the imbibed fluid while only those 
capillaries of Part 2 which abut against larger ones 
are permitted to saturate. Thus, the wetting-fluid 
saturations of Parts 1 and 2 will differ during imbibi- 
tion. The wetting-phase saturation of Part 1 will 
be the same during drainage and imbibition, whereas 
that of Part 2 will be lower during imbibition because 
of the by-passed non-wetting phase. Fig. l(e) com 
pares a model saturation distribution which results 
from drainage with that which resultsfrom imbibition 
to the same capillary pressure. 

The relative permeabilities at any stage of 
the imbibition procees can now be derived from 
the model by determining the area available to flow 
for each of the phases occupying the medium. (The 
expression for ky(jmb) 204 kypiy(imb) are developed 
in the Appendix, and k,,,q,) is shown to be a lower 


limiting value of & y¢imb)-) 


MODEL CHARACTERISTICS 


RELATIONSHIP BETWEEN DRAINAGE AND 
IMBIBITION SATURATION 


The relationship between the drainage and imbibi- 
tion saturation for equal values of k,,,,, which is 
established by the model is given by Eq. 25 *. 





*The development of the formulas expressing the capillary- 
pressure and relative permeability behavior of the imbibition 
model may be found in the Appendix. 
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The family of curves of Fig. 2 shows the imbibition 
saturation vs the drainage saturation with the S,; 
value as a defining parameter. The imbibition model 
is linked to any prototype whose fluid distribution 
is controlled by capillary forces by making the 
model capillary-pressure curve the same as that 
of the prototype. 


MAXIMUM OIL RECOVERY 


The maximum quantity of oil recoverable by ideal 
imbibition conforming to model requirement is equal 
to 50 per cent of the oil initially in place. This 

. . * . 
value is determined by making S,g, = 1 in Eq. 
25. 


DRAINAGE AFTER IMBIBITION 


Let the medium reach a saturation S;,,, through 
an imbibition process. The pores with radius 7, to 
ry are filled by the wetting fluid, except for those 
pores where the non-wetting phase is blocked. 
Assume that at this point the process is reversed 
and that desaturation begins. The biggest pores 
filled with wetting fluid, which are the pores with 
radius r;, and emptied first. At this stage, the fluid 
distribution in all the pores with radii 4,1 to Tp 
has not been modified and is exactly the same as 
it was before the wetting fluid was imbibed into 
pores with radius r;; the permeabilities after desat- 
uration and before imbibition are the same. This 
description of desaturation after imbibition is inde- 
pendent of r; and the relative permeability curves 
which describe desaturation after imbibition are 
the same as those governing the imbibition process. 
The behavior taught by the model conforms to the 
observed stability of relative permeability curves 
after a cycle of oil displacing water and water 
displacing oil.5»6 


RESISTIVITY INDEX 


The resistivity index can be shown to be given 
by 

2 -(l+e€) 

I= Simp acted nae Kine 


where ¢ is a function of S;_,4, smaller than or equal 
toone. The formula for the resistivity index can be 
used to check the consistency of the model. ¢ can 
be computed from experimental measurements of 
I and introduced into the computation of relative 
permeability to the wetting phase. These computed 
values, in turn, can be checked against experimental 
results, 


MODEL APPLICATION 
IMBIBITION RELATIVE PERMEABILITY 


Although the assumptions necessary to build the 
model are too sweeping and arbitrary for the results 
to be rigorously applicable to real reservoir rock, 
the agreement between measured and predicted rel- 
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ative permeability is surprisingly good. Fig. 3 
(a, b and c) shows three sets of measured curves 
from different sources5+7 with the corresponding 
computed values. Here, the agreement is excellent. 
In general, as the rock characteristics shift toward 
those of unconsolidated sands (i.e., toward high 
permeability and low S,,; values), the conformance 
between predicted and observed behavior becomes 
poorer. From the experience acquired in checking 
the model predictions against available data, the 
agreement seems to be generally good for con- 
solidated sands characterized by S,,; values in 
excess of 15 per cent. 

Fig. 2 shows the model definition of the water 
Saturation for an idealized porous medium which 
has undergone imbibition (water invasion) from 
the irreducible water saturation. It is a graphical 


picture of the solution of Eq, 25. The S,,; value 


designating each curve is the parameter used to 
characterize the pore geometry. Hence, if one can 
measure or otherwise estimate S,,;, one can define 
in a few minutes the corresponding imbibition oil or 
gas relative permeability curve. 

The practice of working wholly from the S,; value 
of a sand is acceptable only in the absence of any 
measured flow data. When gas relative permeability 
curves are available it is recommended that they be 
used to define the imbibition k,, curve from Fig. 3. 
Use of measured gas relative permeability data com- 
pensates, at least in part, for the deviation of the 
actual pore-size distribution from the distribution 
which is assumed when working from S,,; values 
alone. 

It is interesting to note that, if the relation be- 
tween S.,*imp and S,*%g, given by Eq. 25 is ac- 
cepted, the equation for the non-wetting drainage 
relative permeability (Eq. 32 of Ref. 3) may be 
transformed into an approximation for the imbibition 
relative permeability. Thus, 


* *2 
m O5-Sy, imb 
keg (imbibition). ~ . O5~ I-Sw imp 


(3) 


RECOVERABLE OIL. 


It is evident that the model value of 0.5 for the 
fractional recovery of oil in place results from the 
assumptions of random pore interconnection and of 
trapped non-wetting phase. These assumptions are 
substantiated by the fact that they lead to predic- 
tions agreeing closely with observations made on 
many consolidated media. 

The development of both drainage and imbibition 
models is based on the premise that the equilibrium 
saturations, at which the effective permeabilities 
become (approach) zero, depend on the pore-size 
distribution. Since the assumed interconnections in 
the models are certainly oversimplified, the behavior 
of a reservoir sand having a particular S,; value 
might depart from that predicted by the model. This 
will be particularly true near the equilibrium satu- 
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rations when the extremes in the pore sizes dominate 
the behavior of the fluid. Thus, as an imbibition or 
drainage process approaches the limiting values, 
one expects and observes differences between 
model and sand behavior. 

It is difficult to assess the agreement between 
the residual oil saturation predicted by the model 
and that which actually occurs. Rough estimates of 
oil recoverable by water flood quite commonly match 
the 50 per cent figure predicted by the model.8 
Nevertheless, firm field or laboratory evidence based 
on data which meet model requirements are hard to 
find. Much of this difficulty stems from the infre- 
quency with which reliable S,,; values are reported. 
The curves of Fig. 3 show the good agreement 
which would obviously result were steady-state 
relative permeability tests used to estimate re- 
siduals. In such tests, the saturation changes are in- 
duced in a gradual stepwise manner so the condition 
for an imbibition saturation distribution are favorable. 
In a flood on a small test plug made at a high 
pressure gradient, the conditions favorable for 
imbibition are doubtful and poorer agreement between 
observed and predicted values is expected. 

The laboratory data of several investigators 9-12 
who used widely divergent displacement techniques 
and the locus of So, predicted by the model are 
plotted on Fig. 4 as S,,;.vs Sg, The obvious scatter 
of the data seems to confirm rather than condemn 
the usefulness of the model prediction. The plotted 
data include tests by Welge 1° and Stahl 9 which 
were made under complete capillary control, flood 
tests on both sands and limestones and flood tests 
on cores containing free-gas saturation induced by 
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FIG. 2~MODEL RELATIONSHIP BETWEEN DRAINAGE 
AND IMBIBITION SATURATIONS AS A FUNC- 
TION OF S,,j. 


both external and internal gas drives. The predicted 
curve is about as good a ‘fit’? as could be expected. 
Welge cautioned that his data could not be repro- 
duced to better than 10 per cent V,. Where flood 
tests were conducted on cores containing free gas, 
the plotted S,, values on Fig. 4 are values of ‘‘one 
minus the final water saturation achieved’’. 

It may be worth emphasizing that the percentage 
of oil recovery predicted by the model is the oil 
recoverable by idealized imbibition into a completely 
water-wet system. Departure of the surface wetting 
from the ideal water-wet condition may change the 
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recoverable oil even from a sand whose pore geom- 
etry conforms to model requirements. Departure of 
the real distribution of pore sizes from the idealized 
distribution of the model will also cause the ob- 
served permeability to approach a zero value at a 
different limiting saturation. 


LIMITS ON THE ROLE OF WETTABILITY 


The water relative permeability for the imbibition 
cycle departs only slightly from that of the wetting 
phase during drainage. The drainage k,, Curve is 
shown in the Appendix to be the lower limit of ky, 
of an imbibition process. Although little work has 
been done as yet on the behavior of the wetting 
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phase, the imbibition k,,, curve derived from the 
model lies very close to the drainage k,,, curve. 
The conformance between measured drainage and 
imbibition k,,, is so good that it seems justifiable 
for present needs to use the much simpler drainage 
curves for the wetting phase. Experimental evidence 
of the similarity between drainage and imbibition 
ky, Curves is shown in Refs. 2, 5 and 7. 

If this suggestion is adopted, it is easy to con- 
struct ky /k» curves for a consolidated sand having 
a known S,,, value. The curves can be constructed 
for both drainage 3.13 and imbibition processes. If 
water invades the sand, its limiting behavior when 
completely water-wet must be characterized by the 
imbibition k,,/k,,, curves of Fig. 5. These are 
labelled water-wet to signify the limiting flow be- 
havior expected during imbibition of water. Were the 
sand completely oil-wet, water invasion would be 
a drainage process and the limiting curves labelled 
“oil-wet’’ would apply. Thus, for a consolidated 
oil sand subjected to a water flood, the curves of 
Fig. 5 define the limits between which flow behavior 
can be influenced by wettability. Similar curves 
for other S,; values readily can be constructed 
to show clearly the increased sensitivity of water- 
oil flow behavior on wettability as the S,,; values 
increase. 

The limiting ratio curves are shown for two sands 
characterized by different S,; values; S,; = 0.30 
and 0.15, respectively. If the sand having an S,,; 
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value of 0.30 were to change from a water-wet 
condition to a completely oil-wet condition, the 
hkyo/kyy tatio and the 50 per cent saturation value 
would decrease to 1/600 of its water-wet value. 
In contrast, a sand characterized by an S,,; value 
of 0.15 would experience a much less-marked change 
in flow behavior due to a like change in wettability, 
At 50% saturation, the k,,,/k,,,, ratio would decrease 
to only 1/30 of the water-wet value. 

A point of interest lies in the interpretation of 
flooding tests which may have been conducted on 
test plugs of intermediate wettability. Oil residuals 
shown by flood tests on plugs which had been made 
pattially oil-wet by mud filtrate or surface handling 
may be erroneously low, as shown by Fig. 5(a). 
However, the same conditions in a test plug of low 


permeability (Fig. 5b) might indicate an erroneously 


high residual oil saturation. A more detailed picture 
of the effect of wettability is shown in Fig. 6. A 
marked difference is evident between the oil relative 
permeability curves measured on oil-base mud cores 
before and after extraction. For comparison, the 
limiting oil-wet and water-wet behavior predicted by 
model theory is also shown. Note how the inter- 
pretation of the intercept, k,, = 0, as a residual oil 
can be misleading when the condition of the rock 
surface is uncertain. 

It is conceivable that the practical use of the 
imbibition model to estimate recoverable oil may 
lie in the limiting ratio curves. By the use of these, 
the effective residual oil saturation would not be 
decided by the ultimate saturation achievable but 
by the permeability ratio beyond which further oil 
recovery becomes unattractive. 


EVIDENCE OF RESERVOIR WETTABILITY 


A search of company files and the industrial 
literature has failed to uncover an example of oil- 
water k,,/k,,, curves which lie near the oil-wet 
limits. The experimental data plotted on Fig. 5 
show some of the correlations. 6,14,15 Although the 
Sig Values of the test sands were rarely known 
exactly, the approximate values assigned are 
certainly close to those for which the correlating 
model curves are shown. Even the Bradford sand, 
so often described as oil-wet, yields ky,/k,,, curves 
very near the water-wet limit. A few data gathered 
by Caudle, et al,14 suggests a sand of appreciable 
intermediate wettability. However, the S,,; values 
reported for this sand was only 1.8 per cent V,, 
an unseemly low value for the other properties. 
Some oil-water relative permeability data for high- 
permeability cores, which were badly contaminated 
with oil-base mud filtrate, fall midway between the 
oil-wet and water-wet limiting curves. 


OIL RECOVERY IN THE PRESENCE OF 
FREE-GAS SATURATION 


If one extends the model behavior to a three- 
phase system under the conditions that water is 
the dominant wetting phase with respect to oil and 
oil with respect to gas, then the fluid distributions 
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resulting from water imbibition will again depend on 
relative pore sizes. The residual saturations can 
be predicted by reasoning similar to that developed 
in the Appendix for two-phase imbibition. 

The maximum water saturation will be 


>05 + oxi a 


The residual oil saturation will be 
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W, imb! max 








where S,; is the initial oil saturation. 
The residual gas saturation fills the remainder 


of the space or 
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The fractional oil recovery in terms of the initial 
gas saturation Sp; is 


Sgi 


No. 
2 1-S\,) 


+ 4 
N 4 


Based on such a model, the group S,;/(1-S,;) 
governs the recovery of oil by water imbibition into 
water-wet sands containing a free-gas saturation. 
On Fig. 7(a) the fractional oil recovered is plotted 
vs the parameter S,;/(1-S,;). Some control of oil 
recovery by this functional group is indicated by 


the Holmgren 11 and Morse data which is compared 


to the curve predicted by Eq. 12. Data from the 


Kyte12 paper show no such correlation on this type 
of plot. 


RESIDUAL GAS SATURATION 


When the model predictions of final saturations 
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are checked against observed residual gas satura- 
tions, the agreement is good only if long times 
for imbibition are prohibited. 

The gas residual after a prolonged time for imbib- 
ition will be considerably lower than that predicted. 
The present interpretation of this apparent discrep- 
ancy is that, if imbibition tests on plugs were con- 
ducted at high pressures (where diffusional effects 
are relatively small), the residual gas saturation 
would have been considerably higher. The paper by 
Gardner, Messmer and Woodside! reported the 
possibility that significant quantities of trapped 
gas may escape by diffusion at low pressures. 
Before a final decision can be reached regarding 
the true meaning of an experimental or computed 
residual gas saturation, this uncertainty must be 


resolved. 

Fig. 7(b) shows residual gas saturations computed 
by Eq. 11 compared with observed values from flood 
tests in the presence of free gas.11,12 Data from 
the Kyte paper 12 are plotted only for tests in which 
mobile gas was initially present. The data from the 
Holmgren paper !! reflect the presence of trapped 
gas, i.e., gas which was already immobilized when 
the water flood started. 


CONCLUSIONS 


A mathematical image(model) has been proposed 
to describe the behavior of a porous medium when 
imbibing a wetting liquid. This model is based on 
the concept that the non-wetting phase is gradually 
blocked during imbibition. 

A theoretical relationship between behavior ex- 
hibited during drainage and imbibition has been 
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obtained from the model. This relationship has 
been applied to show that the relative permeability 
characteristics during imbibition can be determined 
directly from those observed during drainage, or 
from the drainage capillary-pressure curve. 

Application of this model in conjunction with its 
predecessors for the drainage process suggests 
the following. 

1. The degree of oil or water wetness has less 
influence on the flow behavior in sand as its per- 
meability increases, (S,,; decreases). 

2. In lieu of measured oil-water relative permea- 
bility curves, the limiting oil relative permeability 
curve can be determined from a measured or esti- 
mated S,,,; value. The preferred situation would be to 
have a measured keg curve from which to define the 


ky. limits. 

3. A practical limit on the oil recoverable by water 
flood in consolidated rock will be one-half that 
originally present. 

4. The gas volume trapped by water invasion of 
a gas-oil complex may be estimated. 


NOMENCLATURE 


= resistivity index 

= permeability 

= total pore volume 

= pore volume produced 

capillary pressure 

fluid resistivity 

wetting-phase saturation 

irreducible wetting-phase saturation 

non-wetting phase saturation 

= reduced wetting-phase saturation, equal 
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. APPENDIX 
THE COMPUTATION OF FLOW CHARACTERISTICS 
PERMEABILITY 


The method presented here for the computation 
of permeability closely parallels Marshall’s* treat- 
ment of the problem for the drainage process. It 
provides a foundation for the application of the 
model to the imbibition process. 

When the Faces A and B of the capillaric model 
described in the body of the paper are rejoined at 
random, a series of necks of different shapes and 
sizes are created which will govern the flow through 
the medium. The flow through a mean area of these 
necks can be considered as the flow through a 
circular pipe with an equivalent cross section. 
In computing the mean area, the following assump- 
tions are made, 

1. When the fit of one pore to the next is perfect 


(the perimeters do not cross), the cross-sectional 


area of a neck is taken to be that of the smaller 
pore. 


2. The necks are of circular section. 

3. The effective flow path is not increased by 
imperfect alignment of the necks. 

The widening of the flow path within pores is 
neglected and the resulting error is considered to 
be compensated by those introduced by assumptions 
2 and 3. To determine the equivalent area of flow, 
consider the pores as composed of n classes of 
equal total pore area, each class having pores of 
one size only with radii ry > rg >73>... >t. 
If d*a is the total area available to flow, the area 
of each class is d*a/n, When the Faces A and B 
are rejoined at randon, the area occupied by each 
class of pores on Face A will be in contact with 
equal areas of all classes of pores on Face B, 

Let / be a factor expressing the partial abutment 
of a pore class against the solid matrix. Thus, 
when two pores with radii r; and r; abut against 
each other (r; > r,), the area of the neck open to 
flow will be far. The total area of flow through 
necks with area frr2; is composed of: 

1. j areas /(*a/n2) of A with pore radii bigger 
than or equal to r; abutting against pores of B with 
radii rj. 

2. (7 — 1) areas f(d*a/n2) of B with pore radii 
bigger than 7; abutting against pores of A with radii 
Tr 

The total area of flow a7 is then equal to 

a (igo 

f n2 


For the whole medium, the area of flow is obtained 
by summing Eq. 8 from j = 1 to j =n. It can be 
shown that, when a fluid with viscosity p flows 
under a pressure gradient P through x, capillaries 
with radius ry, x2 capillaries, with radii mr, x, 
capillaries with radius r,, such that x; r2,= x72, 
=...=x, 172, the flow per unit area is proportional 
to the arithmetic average of the component areas. 
Then one can write for the flow per average unit 
neck area, 
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The permeability is then given by 
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This formula can be applied to compute the re- 
lative permeability to one phase by considering the 
pores filled with the other phase as part of the 
solid matrix. When the pores 7 to r;.1 are filled 
with non-wetting fluid, ky, is given a 
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If { is taken equal to *S,%, in Eq. llor to 
p*S% in Eq. 12, the formulas become 
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RELATIONSHIP BETWEEN DRAINAGE AND 
IMBIBITION SATURATIONS 

Consider the medium at irreducible saturation 
and let the capillary pressure be decreased to 


pT aliditnimancinnnenitlidibiea (15) 


cn” oF, 


At this pressure, all the pores with radius r, of 
Part 2 abut against pores of Part 1 with radii bigger 
than or equal to r,, and the saturation of Part 2 is 
the same as that of Part 1. 
Sw,27 Sw,imb = Swi * ASw 
where AS, is equal to (1 - S,;) /n : 
At the same time, the saturation of Part 1 will 
be 


Swit . Swidr = Syjp +t AS, 


(17) 


The next step is to decrease Po, to Po, 9.1- 

The area of the class with radius r,., of Part 
2 which will abut against pores of Part 1 with 
radii bigger than or equal to rg_; is[(n—1) /n] AS,,. 
These pores will be filled with wetting fluid, while 
the pores with area (1/n)AS,, which abut against 
pores with radii r, will be by-passed. Then, for 
Part 2, 


S =S), + AS,+ (I-t) AS, 


w,imb 


and for Part 1, 


Swidr = Swi * 20S y 


This process can be continued until pores with 
radius rj are filled. The saturation of Part 2 will be 
S ee 
= Sy + BSy [14 u-4)] 
w,imb wi w ree n ____(20) 
and the saturation of Part 1, 
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Swdr = Swit (n-j) Sy 


When the number n of classes of pores goes to 
infinity, Eq. 20 takes an integral form. Writing j = 


an, o = x and + = dx gives, after integration, 


| - Swi 
Sw, imb = Swi * 2 (I- 


The saturation of Part 1 is computed in the same 
way and is found to be 


Sw,dr = Swi 2 (1=Sa5) (1 - 


Elimination of a and introduction of 


s* 7 Sw, imb ~ Swi Oe 


w,imb 1-Syj 





gives 
* * 
Sw,imb * Sw, dr ~ 


RELATIVE PERMEABILITY (IMBIBITION) 


The mechanics of imbibition will be examined 
hereafter to derive formulas analogous to drainage 
Eqs. 13 and 14 for the imbibition process. 

The permeability to one phase was computed by 
assuming that the pores full of the other phase 
are part of the solid matrix. During the imbibition 
cycle, a certain number of pores of Part 2, full of 


non-wetting fluid, are by-passed and abut against 
pores of Part 1 which are full of wetting fluid. The 


by-passed pores do not contribute to the flow of the 
non-wetting phase, and for the computation of non- 
wetting relative permeability Part 2 can then be 
assumed to have the same saturation as Part l. 
This means that the medium is assumed to be at the 
saturation of Part] or, in other words, that it arrived 
at its present state through a drainage cycle. How- 
ever, its permeability to the non-wetting phase 
corresponds to the true saturation of Part 2. These 
remarks are valid whether the drainage relative 
permeability to the non-wetting phase is computed 
or measured and provide a way of relating imbib- 
ition properties to drainage ones. If the permeability 
is computed, its value is given by Eq. 13. 


g*>s5,%5 " 


| ' 2 
rw kK (2i-Wr yj 3) 
The permeability to the wetting phase is estimated 


as follows. If Part 2 of the slice had the same 
saturation as Part 1, ky, would be given by: 


7 [2i-c2j-n)] 7, 


- fej 





As a certain number of pores of Part 2 are blocked 
and do not participate in the flow, the expression 
for ky, changes. The number of areas is reduced to 
1424+34...4¢(n-j+D=%[(n-j+D(n-j+ 
2)}, and their flow characteristics are distributed 
as follows: 

One area has necks with cross-sectional area 
fur; 

Two areas have necks with cross-sectional areas 
fur uni and 

n —j + 1 areas have necks with cross-sectional 
areas {nr %. 

The value of k,,, is then 


2% % 
k cat is f@ Sy imp 


rw,imb K 





ae oe 
85 (n j+ W(n-j+2) 


y [i-i-1)]r, 


The value of { which was taken equal to #*S%, in 
the drainage process must be re-evaluated. Indeed, 
since the wetting-phase saturation of Part 1 is 
greater than that of Part 2, the probability of having 
a common area of overlap for pores full of wetting 
liquid is greater than if Part 1 had the same satura- 
tion as Part 2. Therefore, { must be written 


where F* is a function of S%, jm, such that 


and ky, imb is given by 


x3 x2_x% 
© FS w,imb 





S(n-jei(n-j+2) 


For n going to infinity, k,,, takes the following 
integral form 
*3 2 * ~2 
"y x F 
Krw, imb = Sw, imb or) 


w,imb 


os . * 
w,imb Sw, imb -S . 
‘ me ee 
p 2 
) c 


The relationship between ky, jmp and ky, gy follows 
from the observation that the permeability to the 
wetting phase, on drainage, is given by 

fe 

Kew, dr = i ‘, dr 

’ K ? 

Sy * 

sia Sw, dr ~$ 


= dS"___(31) 
c 


* 


0 


and that for S% ims = S %,de 


oe 48 


F 

as ( solihesinds (32) 
rw, imb rw,dr lg 
From the definition of F*, the factor F*/S%, igp is 
greater than or equal to one and 

Kew, imb 2 
kw, dr is then lower limit of ky) jmp 

The value of F* cannot be deduced from theo- 
retical arguments only, and experiments are re- 
quired to determine the error made by writing k,, ;, 
= ky de kk 
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ABSTRACT 

This report concerns the microbial flora found 
throughout the surface facilities of six water-injec- 
tion systems in Texas and Oklahoma, Each system 
is described in detail and water quality data are 
presented when available. 

P seudomonads and sulfate-reducers (Desulfovibrio) 
were the predominate organisms found in the systems. 
The genera Sphaerotilus, Bacillus, Achromobacter, 
Micrococcus, Clostridium, Flavobacterium, and 
Sarcina were also isolated in significant numbers. 
Molds, iron bacteria, sulfur bacteria and soil bacteria 
were found less frequently in the systems. 

It is concluded that a complex microbiological 
flora exists in oilfield water-injection systems. A 
tendency for waters to deteriorate in passage through 
solid bed filters is also noted. 


INTRODUCTION 

The injection of water into a subsurface formation 
is often accompanied by numerous problems, the 
origin and nature of which are only partially under- 
stood, During the past few years, investigators have 
indicated that bacteria are involved in some of these 
problems even though very little specific information 
is available concerning this matter. 

It is surprising that investigators working on 
microbiological problems in water injection have 
never made a complete study of the flora of these 
systems. The logical course to pursue is that of 
identifying the microbial flora and then examining 
each organism individually and in combination with 
others to determine if they produce a particular 
problem. 

The objective of the present investigation is to 
identify the microbial flora in water-injection systems. 
In addition, the quality of the.water in most of the 
systems was determined and some observations made 
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on the extent to which passage through surface- 
handling facilities affects water quality. 

The information derived from this study can serve 
as a foundation for additional basic investigations, 


EXPERIMENTAL PROCEDURE 

The microbial flora of six water-injection systems 
were studied at 50 sampling points. These systems 
are experiencing problems which may be attributable 
to micro-organisms, and two were being treated with 
antimicrobial agents at the time of this study. 

The samples were obtained by inoculating 1.0 ml 
of the water onto the different media at the sampling 
point, The media used for primary isolation included 
nutrient agar, Strokes agar, M-10-E medium and semi- 
solid thioglycollate medium. 

The cultures were returned to the laboratory, 
incubated at 30°C until growth occured; then, pure 
cultures were obtained from individual colonies. 
The pure cultures were identified according to 
standard bacteriological procedures. Those organisms 
not fitting identification criteria exactly, even though 
exhibiting close similarities to known species, are 
listed as unidentified organisms. 

Membrane filters were used to make water-quality 
tests at intervals before and after collection of the 
microbiological samples, The filter discs were 47 
mm in diameter with a 0.45-micron particle-size re- 
tention, and the tests were done under closed con- 
ditions at 20 psi. Membrane flow rates are plotted 
as increasing filtrate volumes, and the total solids 
contents are given when available. 

An effort was made to sample at each place in a 
system where a change in environment occurred. 
In this manner, the effect of passage through filters, 
tanks, hay sections, etc., on the microbial flora and 
water quality could be observed. 


RESULTS 

SYSTEM A 
This waterflood system (Fig. 1) is about four years 
old and daily pumps approximately 110,000 bbl of 
oil-free sour water from a 3,700-ft deep reef. The 


7i 
















water passes through two sets of flue gas scrubbers however, there was less variety in the latter system, 
which lower the hydrogen-sulfide content from 250 Pseudomonads and sulfate-reducers were prevalent 
to 20 to 30 ppm. Chlorine is added to remove the in the recycling systems and at the injection wells, 
residual hydrogen sulfide, and this treatment results The membrane data obtained from the supply water 
in the formation of some colloidal sulfur which portion of the system may be found in Fig. 2. The 
cannot be removed efficiently by the anthrafilt water is of good quality through the first scrubber, 
filter. The supply water is treated with an organic however, it deteriorates during passage through the 
corrosion inhibitor. There is no perceptible change second scrubber. The filtered water has a high 
in quality after passage through a maximum of 11 initial flow rate, which falls off rapidly to a value 
miles of distribution lines. below that of Curve C. Further deterioration is 
Chemical analysis of the supply water gave the evident at the injection well. 
following information; pH, 6.85; sodium, 776 ppm; The membrane data (Fig. 3) on the brine recycle 
calcium, 870 ppm; magnesium, 212 ppm; bicarbonate, system shows that the water is of poor quality and 
403 ppm; chioride, 2128 ppm; sulfate, 1350 ppm; and that coalescing, separation and bed filtration are 
hydrogen sulfide, 250 ppm. of very little benefit. Undoubtedly, appreciable 
The system also handles several thousand barrels quantities of insoluble materials are removed, but 
of produced brine per day in two recycling stations. the solid filter cake which develops on the membrane 
The ratio of recycled brine to supply water is cur- disc (and probably at the wellbore) is highly resistant 
rently less than 1:7. The brine is sometimes pro- to the passage of water. The poor quality of water 
duced sour, expecially in the older parts of the in the recycle systems is due primarily to hydro 
flood, but it is generally considered to be sulfide- carbons; however, bacteria and scale deposition 
free in the native state. The injection wells yield may be responsible for the deterioration noted as 
appreciable amounts of iron sulfide when they are the water passes through both systems. 
backwa shed. A core taken from an injection well in this system 
Table 1 shows that the microbial flora in this was cut to 1 in. diameter X 1} in. long, mounted in 
system is quite extensive. In general, the same plastic and placed ina stainless-steel flow cell. 
organisms were found in the oil-contaminated brine The cell was connected to the discharge header so 
recycle systems as in the oil-free supply water, that water flowed through the core under a pressure 
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FIG. 1 — SYSTEM A. 
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TABLE 1 — MICROBIAL FLORA OF SYSTEM A 
Site |-Gas Scrubber No. 1 Inlet Site 9—Coalescer Inlet (Station 1) 








Alginomonas spp. 
Pseudomonas putida 
Achromobacter spp. 
Desulfovibrio desulfuricans 






Site 2~Gas Scrubber No.2 Outlet 


Beggiatoa spp. 

Bacillus mascerans 
Pseudomonas oleovorans 
Desulfovibrio desulfuricans 





Site 3—Filter Outlet 


Vitreoscilla spp. 
Leptothrix spp. 
Pseudomonas boreopolis 














Site 4—Discharge Header 
(Core Outlet) 
Sarcina auriantica 
Proteus spp. 
Achromobacter spp. 













Site 5—Injection Well 








Escherichia freundii 
Pseudomonas spp. 
Desulfovibrio desulfuricans 









Site 6—Injection Well 


Nocardia spp. 
Thiobacillus spp. 
Pseudomonas fluorescens 
Sphaerotilus natans 


Site 7—Injection Well 





Cellulomonas aurogena 
Thiobacillus spp. 
Desulfovibrio desulfuricans 


g Site 8—Injection Well 
- Pseudomonas striata 


Bacillus mascerans 
Pseudomonas boreopolis 





Thiobacillus spp. 
Thiobacillus spp. 
Alginomonas nonfermentans 
Mycoplana dimorpha 
Desulfovibrio desulfuricans 


Site 10—Coalescer Outlet 
Micrococcus flavus 
Achromobacter liquefaciens 
Desulfovibrio desulfuricans 





Site 11—Filter Inlet 
Sarcina lutea 
Micrococcus varians 
Desulfovibrio desulfuricans 





Site 12—Holding Tank Outlet 


Mycoplana dimorpha 
Bacillus cereus 
Pseudomonas ovalis 
Sphaerotilus natans 
Desulfovibrio desulfuricans 





Site 13—Coalescer Inlet (Station2) 





Corynebacterium spp. 
Nocardia opaca 

Proteus rettgeri 

Sarcina auriantica 

Bacillus cereus 
Desulfovibrio desulfuricans 


Site 14—Coalescer Outlet 
Micrococcus conglomeratus 
Pseudomonas ovalis 
Pseudomonas arvilla 
Desulfovibrio desulfuricans 





Site 15—Skimmer Outlet 


Bacillus cereus 
Thiobacillus nevellus 
Vitreoscilla spp. 
Leptothrix spp. 








Desulfovibrio desulfuricans Desulfovibrio desulfuricans 





Site 16—Holding Tank Outlet 





Aerobacter aerogenes 
Escherichia coli 
Thiobacillus spp. 
Pseudomonas aeruginosa 
Pseudomonas fluorescens 
Pseudomonas putida 
Desulfovibrio desulfuricans 








of 700 psi. Based on the Darcy equation, the core 
showed an initial permeability of 96 md, 40 md 
after eight hours, and 37 md after 16 hours, indicat- 
ing that the permeability declined slowly following 
an initial rapid decline. This permeability stability 
at a fairly high figure correlates with the general 
lack of injection problems in portions of this system. 
Injection problems will probably increase as the 
production of brine increases unless more efficient 
oil removal is accomplished. 


SYSTEM B 

This waterflood system (Fig. 4) is two years old 
and daily pumps approximately 90,000 bbl of water 
(temperature, 58° to 60° F) from alluvial beds 


adjacent to the Arkansas River. The water is pumped 
through a coated line to three injection stations, and 
at two of these stations it is mixed with produced 
brine. There is bare steel in the section containing 
produced brine, and black water has been observed 
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FIG. 3 — SYSTEM A — BRINE RECYCLE STATION 2. 


in this section and in the injection wells. The 
system was under treatment with a quaternary amine 
bactericide at the time the microbiological study 
was made. Chemical analysis of the water gave the 
results shown in Table 2. 


Table 3 shows that pseudomonads were the pre- 
dominant organisms in the system. The bactericide 
was doing an excellent job of reducing the variety 
of flora coming from the supply wells since only 
one organism was isolated from the treated water. 

The sulfate-reducers were found primarily in the 
brine portion of the system, but they may be working 
upstream since they were isolated at Site 3. The 
supply wells contain bare steel, and it is possible 
that they may start producing black water in the 
future when the sulfate-reducers reach this point. 

Appreciable quantities of iron oxide have built 
up in the pump bowls and transfer lines in the 
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FIG. 4 — SYSTEM B. 


fresh-water portion of this system. These deposits 
were believed by some to be caused primarily by 
iron bacteria (Gallionella), while others thought 
they were due to the presence of oxygen in the 
system. Gallionella disc counts have been as high 
as 20,000/m! at various points in the system. Oxygen 
levels periodically ran in excess of 1 ppm, and 
considerable effort has been made to reduce air 
entry in the supply wells. 

A study was made to compare Gallionella counts, 
dissolved-oxygen content and suspended iron during 
six-months treatment with an organic sulfur biocide 
and three months without treatment. The oxygen 
level was reduced to less than 0.05 ppm, and the 
total iron available for oxidation was 3 to 4 ppm. 
The results showed that there was no significant 
increase in suspended iron coincident with an 
increase in Gallionella to a point too numerous to 
count. 

Fig. 5 exhibits the membrane data accumulated 
during the biocide treatment at sampling points at 
both ends of the system. The sample taken in January 
was just prior to initiating treatment. Water quality 
improved markedly through June and has continued 
to improve since discontinuing biocide treatment. 

Fig. 5 also gives an analysis of the solids 
retained on the filter disc. The acid-soluble portion 
is represented as Fe (hydroxide and/or sulfide) and 
Ca (carbonate and/or sulfate). That portion insoluble 
in acid but ignitable is reported as ‘‘organic’’. 

The brine is passed through a graphilter filter 
ptior to mixing with fresh water, since it is consid- 
erably poorer in quality than the fresh water. This 
practice results in a water of fair quality being 





TABLE 2—CHEMICAL ANALYSIS OF WATER FOR SYSTEM B 


Brine 


Fresh Water (Average Val.) 


pH 7.20 6.60 

Sodium 238 ppm 23,257 ppm 
Calcium 118 ppm 3,565 ppm 
Magnesium 25 ppm 902 ppm 
Bicarbonate 171 ppm 363 ppm 
Chloride 518 ppm 44,591 ppm 
Sulfate 43 ppm 32 ppm 
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TABLE 3—MICROBIAL FLORA OF SYSTEM B 
Site 1 — Treated Supply Well 
Aerobacter aerogenes 





Site 2 Untreated Supply Well 


Sphaerotilus natans 
Pseudomonas putida 
Pseudomonas fluorescens 
Pseudomonas boreopol is 
Aerobacter aerogenes 





Composite Supply Water 





Desulfovibrio desulfuricans 


Inlet to Fresh-Water Storage Tank, Station | 





Desulfovibrio desulfuricans 
Cytophaga rubra 
Flavobacterium ferrugineum 
Pseudomonas fluorescens 
Pseudomonas putida 


Inlet to Fresh-Water Storage Tank, Station 2 





Sphaerotilus natans 
Pseudomonas aeruginosa 
Pseudomonas boreopolis 


Injection Pump Inlet, Station 3 





Alginomonas nonfermentans 


Produced Brine, Above Filter, Station 2 
Desulfovibrio desulfuricans 
Sphaerotilus natans 

Pseudomonas fluorescens 








furnished to the injection wells. Fig. 6 shows that 
the filter is doing a good job of total solids removal, 
but membrane plugging persists. This suggests that 
only the large particles are removed. 


SYSTEMC 

This pressure-maintenance system (Fig. 7) is 
about four years old and daily injects 250 to 1,000 
bbl of brine produced from wells about 7,300-ft 
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FIG. 6 — SYSTEM B — PRODUCED BRINE. 


deep. The temperature of the brine varies from 75° 
to 100°F, and it contains a considerable amount of 
entrained oil in the form of an inverse emulsion. 
The injection well has plugged severely and does 
not respond satisfactorily to acidification. A blend 
of corrosion inhibitor, scale preventive and surfactant 
is added to the water. Chemical analysis of the 
composite water gave the following information: pH, 
6.2; sodium, 20,423 ppm; calcium, 1,120 ppm; 
magnesium, 249 ppm; bicarbonate, 268 ppm; chloride, 
34,042 ppm; and sulfate, 24 ppm. 












































Table 4 shows that this system is extensively con- 
taminated with bacteria and that most of the organisms 
were encountered around the skimmer tank and injec- 
tion well. This system, like the others, contains 
pseudomonads and sulfate-reducers but, unlike the 
other systems, also contains a number of molds. 

The membrane data (Fig. 8) show that the ultimate 
plugging time of the water was not materially improved 
by passage through the conditioner and that an appreci- 
able quantity of nonferrous acid-soluble material 
(mostly calcium carbonate) is present. The organic 
matter does not include oil, which was found to be 
20-ppm above and 10-ppm below the conditioner. 

A core flow test was made at sample Point 6(Fig. 
7) with a l-in. plug of Berea sandstone at 100 psi. 
The initial permeability was 64 md, which was reduced 
to 4.2 md after one hour, and 1.9 md after five hours. 
The permeability of the core was increased to 8 md 
after treatment with acid and organic solvents with 
and without surfactants. 

The data suggest that bacteria and suspended solids 
are responsible for the plugging problems encountered 
in this field. 


SYSTEM D 


This waterflood system (Fig. 9) has been in oper- 
ation for about eight years and daily injects approx- 
imately 3,500 bbl of produced brine and 6,500 bbl of 
fresh water. The temperature of the produced brine 
varies from 75° to 100°F, and the fresh water from 50° 
to 60°F. An aldehyde bactericide is added to the water 
because the system has developed several problems 
such as plugged wells, black water and corroded 
meters. Table 5 shows the chemical analysis of the 
waters. 


Table 6 shows that the system is contaminated with 
Sphaerotilus and sulfate-reducers, which may account 
for the black water. Only one organism was isolated 
from the fresh water, and the variety of flora is not 
great throughout the system. 

Membrane data (Fig. 10) show that mixing good- 
quality fresh water with the brine produces a composite 
water with rapid-plugging tendencies. The standard 
graphilter-pack filter(minimum particle size 1/16 in.) 














HOLDING Tan} —€yp—— © 





ONE HALF 
MILE 











OTHER 
LEASES 


FIG. 7 — SYSTEM C. 
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TABLE 4 — MICROBIAL FLORA OF SYSTEM C 
Site 1—Clear Water 4-Ft Off Bottom of Gunbarrel 

Pseudomonas aeruginosa 

Bacillus megaterium 


Corynebacterium poinsettiae 
Unidentified organism 





Site 2—Sediment 2 to 3 In. Off Bottom of Gunbarrel 


Micrococcus agilis 
Sphaerotilus natans 
Pseudomonas putida 
Pseudomonas oleovorans 
Bacillus sphaericus 
Desulfovibrio desulfuricans 
Two unidentified organisms 





Site 3—Centrifugal-Pump Discharge Between Skimmer Tank and 
Conditioner 





Achromobacter stenohalis 
Pseudomonas oleovorans 
Pseudomonas putida 
Pseudomonas fluorescens 
Mycobacterium marinum 
Micrococcus halodenitrificans 
Penicillium notatum 
Penicillium spp. 

Mucor spp. 

Asperfillis spp. 

Monilia spp. 

Rhizopus spp. 

Desulfovibrio desulfuricans 
Three unidentified organisms 


Site 4—Line From Bottom of Filter 
Bacillus cereus 
Desulfovibrio desulfuricans 
Two unidentified organisms 





Site 5—Discharge Line From Bottom of Horizontal Holding Tank 





Bacillus cereus 

Bacillus cereus, var. mycoides 
Desulfovibrio desulfuricans 
Two unidentified organisms 


Site 6—Line Between Horizontal Holding Tank and Injection Well 





Micrococcus varians 
Micrococcus flavus 

Sarcina flava 

Sarcina lutea 

P seudomonas fluorescens 
Pseudomonas putida 
Pseudomonas oieovorans 
Alcaligenes metalcaligenes 
Cytophaga rubra 
Desulfovibrio desulfuricans 
Two unidentified organisms 





produces an unusually good water which undergoes 
extreme deterioration by the time it reaches the injec- 
tion well. The problems encountered in this system 
may be of microbiological origin since entrained oil 
does not appear to be the major source of trouble. 


SYSTEM E 


This pressure-maintenance system (Fig. 11) has 
been in operation for about five years and daily in- 
jects between 30,000 to 60,000 bbl of sour water 
obtained from the same reef as in System A. A semi- 
polar organic corrosion inhibitor is added (25 ppm) 
continuously at the supply wells, and the water is 
then pumped through an 18-in. bare-steel line to the 
injection plant where it is filtered through 12 parallel 
anthrafilt filters (0.6 to 0.8 mm). Chemical analysis 
of the water gave the same results as in System A. 
It was also noted that the sulfate and sulfide content 
of the water remains unchanged throughout the system. 
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FIG. 9 — SYSTEM D. 


Table 7 shows that Desulfovibrio, Sphaerotilus and 
Flavobacterium species were prevalent throughout 
the system. Bacteriological analysis four months 
prior to this study also indicated that the system was 
extensively contaminated with sulfate-reducers. 

Membrane filter analysis shows that the filters 
(Fig. 12) are effectively removing solids which con- 
sist predominantly of iron sulfide. The water un- 
doubtedly deteriorates before it reaches the injection 





TABLE 5—CHEMICAL ANALYSIS OF SYSTEM D WATERS 


Composite 
Fresh Water __Brine 


pH 7.4 6.8 


15,876 ppm 
2,555 ppm 
604 ppm 
220 ppm 
30,645 ppm 
10 ppm 


191 ppm 
39 ppm 
19 ppm 

232 ppm 
11 ppm 
11 ppm 


Sodium 
Calcium 
Magnesium 
Bicarbonate 
Chloride 
Sulfate 
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FIG. 10 — SYSTEM D. 


wells, but no injection problems have been encoun- 
tered in the porous dolomitic reservoir. 


SYSTEM F 


This salt-water-disposal system (Fig. 13) has 
been in operation for 10 years and daily injects 
approximately 3,000 bbl of produced brine obtained 
from wells that are about 7,000 ft deep. Several 
years ago it was suspected that entrained oil was 
plugging the injection well. Temporary benefit was 
obtained by acidizing the well, and treatment with a 
nonionic surfactant was helpful for several months 
but eventually ceased to be effective. A kerosene 
wash was installed at that time, which increased 
intervals between backwashing and/or acidizing 
considerably. The well is having no injectivity 
injectivity difficulties at the current rate of water 
input, and for this reason the system is not under 
chemical treatment. Chemical analysis of the water 
gave the following results: pH, 7.1; Calcium, 1,935 
ppm; Magnesium, 165 ppm; Barium, 108 ppm; Stron- 
tium, 91 ppm; Bicarbonate, 188 ppm; Chloride, 
55,000 ppm; and Sulfate, 3 ppm. 

The system is not experiencing any trouble that 
can be attributed to bacteria even though Table 8 
shows that it is heavily contaminated with micro- 
organisms. Desulfovibrio species were found in 
every sample taken from the system, but black water 
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TABLE 6 — MICROBIAL FLORA OF SYSTEM D 
Site 1—Producing Well Untreated Site 6—Composite Water Below 
Filter 
Desulfovibrio desulfuricans 
Sphaerotilus natans 
Arthrobacter spp. 








Micrococcus lutea 
Bacillus megaterium 
Arthrobacter ureafaciens 
Arthrobacter citreus 
Bacillus cereus 

Site 7—Composite Water at Clear 

Site 2~Producing Well, Treated Tank 

Sphaerotilus natans 
Sarcina ureae 
Micrococcus flavus 
Unidentified yeast 





Desulfovibrio desulfuricans 
Sphaerotilus natans 
Proteus vulgaris 

Proteus rettgeri 


Site 8—Composite Water, Meter 
Block 
Desulfovibrio desulfuricans 
Alginomonas nonfermentans 
Achromobacter liquefaciens 


Site 3~Tank Battery No.1 
Desulfovibrio desulfuricans 
Sphaerotilus natans 
Unidentified yeast 








Site 4—Fresh-Water Supply Well 
Achromobacter xerosis 





Site 9~Backflow at Injection Well 





Desulfovibrio desulfuricans 
Achromobacter liquefaciens 
Clostridium histolyticum 


Site 5—Composite Water Before 
Filter 
Desulfovibrio desulfuricans 
Sphaerotilus natans 








has never been a problem. The very low sulfate 
(SO,4) content of the brine may account for this 
observation. 

Curve A of Fig. 14 shows that one of the lease 
gunbarrels produces a fair-quality brine. The com- 
posite brine is less filterable, but it is improved by 
the kerosene-wash treatment. Passage through the 
open pit and sand filter results in a very poor-quality 
water going to the injection well. 


DISCUSSION 


Table 9 shows that Desulfovibrio, Pseudomonas, 
Sphaerotilus and Bacillus were the most predominant 
genera found in the systems. These organisms have 
been found previously in injection systems, and they 
have been implicated in hydrogen-sulfide production, 
metal corrosion, slime formation, hydrocarbon oxida- 
tion and plugging problems,'~* 

It is interesting to note that pseudomonads and 
sulfate-reducers (Desulfovibrio) were isolated from 
all systems. These organisms are often found to- 
gether where microbiological problems are encoun- 
tered in the petroleum industry. They have been 
repeatedly observed in waterflood systems, in the 
soil surrounding pipelines and in deteriorated petro- 
leum products.°°? 

Several investigators have suggested that there is 
an interrelationship between pseudomonads and 
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FIG. 11 — SYSTEM E. 
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FIG. 12 — SYSTEM E— SUPPLY WATER. 


sulfate-reducers §»9, The pseudomonads may oxidize 
the hydrocarbons and feed the sulfate-reducers, and 
the latter organisms may protect the pseudomonads 
from antimicrobial agents. 

There is no relationship between the age of the 
system and the type of flora encountered. The pre- 
dominant flora of the newer systems is not significi- 
cantly different from that of the older systems, and 
the number of different species isolated per sample 
cannot be related to age of the system. 

No specific differences in the microbial flora of 
treated and untreated systems could be detected 
Fewer species (per sample) were isolated from the 
systems under treatment with antimicrobial agents, 
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TABLE 7 — MICROBIAL FLORA OF SYSTEM E 
Site 1— Supply Wellhead (Treated Water) 


Desulfovibrio desulfuricans 
Sphaerotilus natans 
Flavobacterium rhenanum 


Unidentified organism 





Site 2— Supply Wellhead (Untreated Water) 
Desulfovibrio desulfuricans 
Sphaerotilus natans 
Achromobacter liquefaciens 
Two unidentified organisms 





Site 3— Composite Water, 18-In. Line Near Scraper Inlet 





Desulfovibrio desulfuricans 
Flavobacterium lutescens 
Micrococcus rubens 
Scopulariopsis brevicaulis 


Site 4— Composite Water Before Entering Filter Pump Station 





Desulfovibrio desulfuricans 
Sphaerotilus natans 
Flavobacterium arborescens 
Unidentified organism 


Site 5— Composite Water Below Filter Station 





Desulfovibrio desulfuricans 
Flavobacterium arborescens 
Pseudomonas aeruginosa 
Bacillus coagulans 


Site 6— Composite Water, Southeast Section of Field Near Injection 
Well 
Desulfovibrio desulfuricans 
Sphaerotilus natans 
Unidentified organism 








and one system receiving a corrosion inhibitor 
contains the greatest variety of flora. 

The nature of the injection water may influence 
the variety of organisms isolated from the systems. 
The greatest number of different species (per sample) 
were found in two systems handling produced brine 
containing appreciable quantities of entrained oil. 

The membrane and solids data are of limited value 
if they cannot be related to the ability of the forma- 
tion in question to take water. A water is of poor 
quality if it plugs a 0.45-micron disc to a rate of 
less than 2 ml/sec; at 1,000-ml throughput at 20 psi, 
however, it may not produce plugging problems. This 
depends on the permeability and porosity of the for- 
mation and the surface pressure. More work of the 
type reported on cores in Systems A and C is needed. 
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FIG. 13 — SYSTEM F. 
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With appropriate core flow data, it should be possible 
to predict much more accurately the anticipated 
injectivity into a formation of known capacity. 
The problem of effective surface water handling 
continues to be largely unsolved in the systems 
studied. The membrane filter tests show that bed- 
type filters generally do not improve water quality. 
These filters remove considerable solids and prevent 
occasional slugs of oil from entering the formation, 
but they do not remove fine solids. In addition, the 





TABLE 8 — MICROBIAL FLORA OF SYSTEM F 


Site 1]—Gunbarre!l Produced 
Water, Bottom Level 
Pseudomonas putida 
Bacillus sphaericus 
Micrococcus lutea 
Desulfovibriodesulfuricans Site 5—Open-Pit Outlet 
seudomonas fluorescens 
Pseudomonas putida 
Arthrobacter pascens 
Bacillus cereus 
Desulfovibrio desulfuricans 


Site 4—Kerosene-Wash Outlet 





Clostridium tertium 
Pseudomonas fluorescens 
Desulfovibrio desulfuricans 





Site 2—Produced Water At Oil- 
Skimmer Outlet 


Pseudomonas aeruginosa 
Arthrobacter simplex 
Arthrobacter ureafaciens 
Sarcina lutea 
Achromobacter iophagus 
Flavobacterium spp. 
Escherichia coli 
Desulfovibrio desulfuricans 





Site 6—Injection Well Below 
Filter 

Pseudomonas aeruginosa 
Pseudomonas boreopolis 
Bacillus megaterium 
Bacillus cereus 
Bacillus polymyxa 
Escherichia intermedia 
Proteus rettgeri 
Flavobacterium spp. 
Arthrobacter spp. 
Desulfovibrio desulfuricans 





Site 3—Kerosene-Wash Inlet 





Clostridium histolyticum 
Proteus inconstans 
Pseudomonas fluorescens 
Bacillus cereus 
Desulfovibrio desulfuricans 





TABLE 9 — SUMMARY OF THE 
MICROBIAL FLORA OF ALL SYSTEMS 
Per Cent Times Isolated From System 
Genus A B C D E_ F_ Average 
Desulfovibrio spp. 42 83 66 100 100 78.6 
Pseudomonas spp. 57 66 16 100 51.0 
Sphoerotilus spp. 42 16 66 0 31.8 
Bacillus spp. 0 66 1% 66 31.6 
Clostridium spp. - - 33 22.0 
Micrococcus spp. 50 1% 1% 20.3 
Flavobacterium spp. 0 66 833 18.8 
Achromobacter spp. 16 1%6 16 16.5 
Arthrobacter spp. 0 0 = 650 12.0 
Sarcina spp. 16 0 16 10.1 
Proteus spp. 0 33 9.3 
Escherichia spp. 0 33 7.5 
Alginomonas spp. W 6.1 
Aerobacter spp. 5.6 
Yeast 5.5 
Thiobacillus spp. 5.1 
Cytophaga spp, 5.0 
Scopulariopsis spp. 2.6 
Penicillium spp. 2.6 
2.6 
2.6 
2.6 
2.6 
2.6 
2.6 
2.6 





Alcaligenes spp. 
Aspergillis spp. 
Corynebacterium spp. 
Monilia spp. 
Mucor spp. 
Mycobacterium spp. 
Rhizopus spp. 
Leptothrix or 
Toxothrix spp. 
Mycoplana spp. 
Nocardia spp. 
Vitreoscilla spp. 
Beggiatoa spp. 
Cellulomonas spp. 
Corynebacterium spp. 


coooooooooowococno 


2.0 
2.0 
2.0 
2.0 
1.0 
1.0 
1.0 
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filters and other surface vessels are excellent 
environments for microbial growth. Diatomaceous 
earth filters remove solids more efficiently, but they 
usually represent a high initial investment and have 
a much lower flow rate per square foot of filter surface. 
Bare-screen type filters or strainers have applica- 
tion in installations where there are only traces of 
hydrocarbons. 

There is a definite need for improvement in oilfield 
water-filtration equipment, but it may not be forth- 
coming until more definite water-quality standards 
are established by the operators. 


CONCLUSIONS 


1. A wide variety of microbial flora is present in 
the six water-injection systems. Four genera con- 
stitute the major portion of the flora; however, other 
previously unreported organisms are also present. 

2. Only one system appeared to be experiencing 
problems that may be due to micro-organisms. 

3. Membrane data are of comparative interest but 
do not relate directly to well injectivity unless 
accompanied by core flow data on the formation in 
question. 

4. Membrane filter tests on the systems show that 
bed filters do not usually improve water quality. 
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Injection Rates—The Effect of Mobility Ratio, Area 


Swept, and Pattern 
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MEMBER AIME 


ABSTRACT 


A method is presented for calculating approximate 
injection rates in secondary recovery operations. 
The method can be applied to cases of unequal fluid 
mobilities, irregular well patterns and boundary 
patterns. 

The steady-state pressure distributions for the 
four flood patterns reported by Muskat and for five 
additional patterns reveal that most of the difference 
in pressure between the injection and producing 
wells occurs in regions around the wells which can 
adequately be described as regions of radial flow. 
This leads to a method of calculating injectivity by 
approximating the flood pattern with radial flow 
elements (or a combination of radial- and linear- flow 
elements for some patterns such as the direct line 
drive), Irregular and boundary patterns can also be 
approximated by radial and linear elements. 

Each of these elements can be described by radial- 
and linear-flow equations and the results combined 
as series flow resistances to give an approximate 
equation for the initial injection rate. The mobility 
ratio does not affect the initial rate; therefore, if 
the well pattern is one of those regular patterns for 
which theoretical rate equations have been derived 
for unity mobility ratio, the approximate initial rate 
equation can be improved by adjusting it to match 
the theoretical equation. The available theoretical 
rate equations are listed, including five new cases. 

As the flood progresses, the injectivity changes 
because in general the flood front will divide the 
pattern into areas of different fluid mobilities. Sim- 
ple shapes can be assumed for the flood front so 
that both areas can be divided into radial- and 
linear-flow elements. 

Radial- and linear-flow equations are applied to 
these elements to account for the change in flow 
resistance behind the front. 

To calculate injection rates after breakthrough, 
it is necessary to know the sweep efficiency at 
breakthrough. Areal sweep data are available in the 
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literature for a number of patterns, and sufficiently 
accurate breakthrough sweep efficiency can be esti- 
mated from these data if it is not otherwise avail- 
able. 

Again, simple shapes for the flood front can be 
assumed after breakthrough, and injection rates can 
be calculated for the remainder of the flood. When 
areal sweep data are available, these data can be 
used to check the calculated injection rates after 
breakthrough. 

Injection rates calculated by this approximate 
method compare favorably with available model data. 


INTRODUCTION 


When the fluid mobilities in the swept and unswept 
regions are equal, the injectivity will not change as 
the flood front advances; and for regular patterns it 
can be calculated by mathematical formulas.! 

When the fluid mobilities in the swept and unswept 
regions are not equal, the injectivity will increase 
or decrease as the flood front advances. In this 
case, the injectivity has not been calculated by 
analytical means for any practical well pattern; and, 
furthermore, scale model> and analog**® results 
have been published only for the five-spot pattern. 

The main object of this paper is to present an 
approximate method for calculating injectivity for 
the case of unequal mobilities. The method can be 
applied to regular, irregular and boundary patterns. 

Before the approximate method is discussed, the 
analytical formulas for equal mobilities will be 
summarized, and the analytical solution for radial 
flow with unequal mobilities will be used to show 
how the injectivity changes as the flood progresses. 


EQUAL MOBILITIES, REGULAR PATTERNS 


In a flooding operation the injectivity, the rate at 
which fluid can be injected per unit difference in 
pressure between injection and producing wells, 
depends on basic physical properties and, in addi- 
tion, on three variables to be treated here—the area 
of the swept region, the fluid mobilities in the swept 
and unswept regions, and the well geometry (pattern, 
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spacing and wellbore radii). 

When the fluid mobilities in the swept and unswept 
regions are the same, the injectivity of a given well 
atray is independent of the first variable, the size 
and shape of the area swept, and is directly pro- 
portional to the single fluid mobility involved. The 
calculation of the injectivity thus reduces to the 
geometrical problem of treating the particular well 
pattern, well spacing and well radii. 

For regular patterns, Muskat! has offered a method 
for finding analytical mathematical solutions and 
has worked out the formulas for special cases — 
the direct line drive, the staggered line drive, the 
five-spot and the seven-spot. Five more cases are 
presented here--the nine-spot with arbitrary corner- 
well to side-well producing ratios, a direct-line- 
drive boundary pattern, two five-spot boundary pat- 
terns and a nine-spot boundary pattern. These nine 
cases are listed in Table 1, which is thus a sum- 
mary of the available results on injectivity of regular 
patterns with unity mobility ratio. The derivations 
of the five new cases are omitted because they 
follow the method of Muskat. 

The five new cases can be used in the develop- 
ment of an approximate method for handling irregular 
patterns. The nine-spot will be used to show the 
extent of approximately radial flow even for large 
differences in corner-well and side-well rates. The 
four boundary cases can be used to aid in estimat- 
ing the angle open to flow in irregular - boundary 
patterns. 


RADIAL SYSTEM, UNEQUAL MOBILITIES 


In each of the cases appearing in Table 1, the 





TABLE 1A — INJECTIVITIES FOR REGULAR PATTERNS 
WITH UNITY MOBILITY RATIO. 
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injectivity is independent of the area swept and 
depends on mobility only to the extent of simple 
proportionality. The only important variable is the 





TABLE 1B — INJECTIVITIES FOR REGULAR PATTERNS 
WITH UNITY MOBILITY RATIO. 
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TABLE IC — INJECTIVITIES FOR REGULAR PATTERNS 
WITH UNITY MOBILITY RATIO. 
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producing well q'. 


(Wells of outer row are all producing wells.) 
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TABLE 1D — INJECTIVITIES FOR REGULAR PATTERNS 
WITH UNITY MOBILITY RATIO. 
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system geometry and, for these restricted conditions, 
analytical solutions are possible. 

In this section another analytical solution is 
considered, this time made possible by geometrical 
simplicity. This solution illustrates the manner in 
which the other variables enter the problem, show- 
ing explicitly the variation of injectivity with posi- 
tion of the flood front for any mobility ratio. 

Consider a radial system with a central injection 
well of radius r,, and imagine that fluid is produced 
uniformly from each point on a circle of radius re. 
Purely radial flow will result, and the front between 
injected and original fluids will be a circle whose 
radius will be called r,. Simple radial-flow equations 


will apply in both the swept and unswept regions, 
In the unswept region, the production rate q will be 
given by 

_ 0.003076 Ay (P¢ - PF ) 


log (%/r) 





Cae es 2 ee ee Se ee (1) 
where 4 is the thickness (in feet), A, is the fluid 
mobility in the unswept region (millidarcies/centi- 
poise), and pe and py are the pressures at the outer 
boundary and at the front (psi). Similarly, in the 
swept region the injection rate i will be 


_ 0.003076 AAs (pF - Py) 
j ,--.res B/D 
log(%/r,) 





ee ee ee ee ae (2) 
where A, is the mobility in the swept region and 
Py is the pressure at the injection well. 

For steady-state incompressible flow, the injec- 
tion rate is equal to the production rate (i = q), and 
these two equations can be combined to eliminate 
py and give an expression for the injection rate in 
terms of the position of the front r¢ which is valid 
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for any mobility ratio. 
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where Ap =p, — py and M=A,/A,. 

In most applications, it is convenient to express 
the variation in injectivity with the progress of the 
front in terms of the initial injectivity. In this case, 
that is, a relative injectivity I, is defined as the 
ratio of the injectivity at any time (given by Eq. 3) 
to the initial injectivity (given by Eq. 3 with rz = 
Task 
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The relative injectivity starts at unity when Te =Ty 
and ends at M when 77 = re (corresponding to com- 
plete sweeping of the area and complete change of 
fluid mobility from A, to As). 

Curves calculated from this equation are plotted 
in Fig. 1 for mobility ratios of 5, 1 and 1/2, and 
for radius ratios ro/r,,, of 160 and 500. These ratios 
also appear in later comparisons of approximate 
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calculations and model results. The relative injec- 
tivity is plotted against fraction of area swept Es, 
rather than against the position of the front. In this 


. - ea 
case, the relation is "f /te = Eg, for ry Kt. 


The relative injectivity changes very rapidly 
during the early part of the flood, and then shows 
a gradual change during the remainder of the flood. 
This early, rapid change followed by a long period 
of slow change is typical of radial systems. In 
complex systems, the rapid change occurs while 
the flood front is moving out from the injection well 
in a near-radial manner. 

In this simple example of complete radial flow, 
the period of gradual change extends to the end of 
the operation because breakthrough of the injected 
fluid into the producing circle occurs at complete 
sweep of the area. In a more complex system, with 
producing wells rather than the producing circle, a 
nose of injected fluid may reach a producing well 
before complete sweep of the area. After this break- 
through point, there again will be a more rapid 
change in the injectivity. 


APPROXIMATE METHOD 


REGULAR PATTERNS 

Analytical equations of steady-state pressure 
distributions in complex flood patterns have been 
derived by Muskat! for the case of unity mobility 
ratio. An examination of these pressure distributions 
shows that most of the pressure difference between 
injection wells and producing wells occurs in areas 
around the wells which can adequately be described 
as regions of radial flow. In flood patterns as com- 
plex as the nine-spot pattern (see the Appendix), 
these regions of near-radial flow still occur around 
the injection and producing wells, 

Porous, fluid-flow model studies? show that, 
even when the mobilities of the swept and unswept 
regions are different, the flood front moves out 
from the injection well in a radial manner during 
the early part of the flood. 

These properties of flood systems suggest that 
the injectivity of any system can be simply and 
quickly approximated by dividing the pattern into 
segments of radial and linear flow and computing 


the total flow by combining the radial- and linear- 
flow equations which apply to individual segments. 

In estimating injectivity by this approximate 
method, the volume of the pattern must be divided 
among the producing wells of the pattern. The volume 
assigned to each of the wells is then approximated 
by regions of radial and linear flow. The injectivity 
can be calculated by assuming that the flood front 
advances radially in the radial regions and linearly 
in the linear regions until breakthrough is ap- 
proached. In this manner, injectivity prior to break- 
through of the swept region into a producing well 
can be estimated without the need of areal-sweepout 
data for the flood pattern. 

As breakthrough is approached, a nose begins to 
grow at the flood front and streaks to the producing 
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well during a short period of time when only a 
small part of the total-pattern volume is swept. As 
further production occurs, the finger opens up and 
the rate of flow from the swept region is nearly 


proportional to the angle this finger makes at the 
producing well. Prats, et al,6 showed that it is the 
existing angle open to flow which controls the 
injectivity. For a given angle, the exact shape of 
the flood front after breakthrough has little effect 
on injectivity. If the flood front after breakthrough 
is approximated by a circle with a sector whose 
angle is proportional to the rate of flow from the 
swept region, the injectivity after breakthrough can 
be estimated from the same radial- and linear-flow 
segments used before breakthrough. 

To illustrate this method of estimating injectiv- 
ity, consider the five-spot pattern shown in Table 
1. The symmetry element of the five-spot can be 
divided into two radial segments as shown in Fig. 2, 

Applying the radial-flow equation to each of the 


two radial segments and combining gives for the 
initial injection rate 


0.003076 ha, ( Dwi ~Pwp) 





/ ,» res B/D 
™ lei ‘a ‘ep 
ene @) os 
9 wi g ‘wp 
ee ee ee oe ees . (4) 
where h =net pay thickness, ft, 


fe; =tadius of the radial segment around 
the injection well, 
tep =tadius of the radial segment around 
the producing well, 
Twi = effective wellbore radius of the injec- 
tion well, 
Twp = effective wellbore radius of the pro- 
ducing well, 
P wi = pressure at injection well, psi, 
Pwp = pressure at producing well, psi, 
A, = mobility of the reservoir fluid, k,/py, 
md/cp, and 
d=distance between injection and pro- 
ducing wells (see Fig. 2). 
When the effective wellbore radii are equal and 
Tei =Tep Eq. 4 becomes 
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which agrees very closely with the exact equation 
for the five-spot pattern given in Table 1 (the only 
difference is between the constants 0.249 and 0.269). 
A still closer agreement could be obtained by ap- 


proximating the pattern area with a lineareflow seg- 
ment between two radial-flow segments. 

On the assumption that the flood front moves out 
in a true radial manner, the injectivity while the 
front is between the injection well and the radial 
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‘oe, 
+ log ya 
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res B/D-psi 
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where M = mobility ratio A,/A,, and the relation be- 
tween the radius of the flood front rj, and the fraction 
of total area covered by the swept region Eg, is 
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When the flood front is between the boundary re 
and the producing well prior to breakthrough, the 
injectivity is 
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FIG. 3 ~ASSUMED SHAPE OF FLOOD FRONT 
AFTER BREAKTHROUGH OF THE 
SWEPT REGION INTO A PRODUC- 
ING WELL OF A FIVE-SPOT PAT- 

TERN. 
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where /¢ = 7 a | peer ee S 


After breakthrough, the flood front will be assumed 
to have the shape shown in Fig. 3. 


The flow rate from each region is proportional to 
the mobility and to the angle open to flow from that 
region. This leads to the relation 


a oe = 
On £+M(i-f) * 
where / is the fraction of the producing stream from 
the swept region, and a is the angle (in radians) 
open to flow from the swept region. 
With values of a from Eq. 10 and data of areal 
sweep vs fractional flow from the swept region, the 
injectivity after breakthrough is 
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When areal sweep vs fractional flow from the 
swept region is net known, an approximation of 
injectivity after breakthrough can be obtained by 
assuming the radius rj to be constant and equal to 
the value obtained from Eq. 9 or Eq. 12 using the 
breakthrough value of areal sweep E,. With a fixed 
radius rj, the angle a determined by Eq. 10 also 
establishes areal sweep as a function of fractional 
flow from the swept region. Patterns for which areal 
sweep data are available show no significant dif- 
ference between the injectivity calculated using 
areal sweep data and the injectivity calculated by 
fixing the radius rj at the breakthrough value of 
areal sweep. 


The approximate relative injectivities for the 
five-spot pattern calculated from Eqs. 4, 6, 8, 10 
and 11 are shown in Figs. 4 and 5 for comparison 
with experimental data reported by Aronofsky, et al,* 
and by Caudle, et al.5 The agreement is good except 
for very low mobility ratios where the effects of 
by-passing, fingering, and saturation gradient across 
the front become more pronounced. For these cal- 
culations, the value of re; was taken equal to rep 
for mobility ratios for which breakthrough occurs at 
areal sweeps of more than 50 per cent of the total 
area. For mobility ratios in which breakthrough 
occurred prior to 50 per cent areal sweep, the value 
of re; was calculated from the per cent areal sweep 
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= 
at breakthrough as re; =d\/ af , and the value of 
2d? 


Tep was calculated as rep = \/ —— —Tei* 


IRREGULAR WELL SPACING OR 
UNEQUAL PRODUCING RATES 


The injectivity of a flood pattern with irregular 
well spacing or unequal producing rates can be 
estimated in a manner similar to that described for 
the five-spot by distributing the drainage area of 
the pattern among the producing wells. The area 
assigned to each of the wells is then approximated 
by radial- and linear-flow segments, and the progress 
of the front is traced through these segments. 

Dividing the area among the producing wells first 
involves dividing the total area into isolated or 
separate elements (corresponding to the division of 
the five-spot into symmetrical elements) and, then, 
dividing this isolated element among the producing 
wells in the element. 

One method of dividing the total area into isolated 
elements for separate treatment is to use as divid- 
ing lines the lines joining all the producing wells 
which are supplied by one particular injection well. 
Fig. 6 shows an isolated element obtained in this 
way. The centrally located injection well supplies 
four producing wells. Lines are drawn connecting 
each of these producing wells with its two neighbors, 
and these lines form an irregular polygon which is 
taken as an isolated element of the total flood area. 

To distribute this isolated element among the 
producing wells, a point is located on the line 
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FIG. 4—COMPARISON OF APPROXIMATE 
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MENTAL INJECTIVITY OF 
ARONOFSKY FOR FIVE-SPOT 
PRODUCING PATTERNS. 


joining each pair of neighboring wells. This point 
divides the line between the wells in proportion to 
the total producing rates of the wells. In Fig. 6, 
for example, suppose that the total well rates have 
the following ratios: q1/q2 = 1, 92/93 = 5/4, 93/94 
= 8/9, and q4/q, = 9/10. Then four points are lo- 
cated as shown in Fig. 6, dividing the four sides of 
the irregular polygon into these ratios. Points are 
located at one-half the distance from Well 1 toward 
Well 2, five-ninths the distance from Well 2 toward 
Well 3, and so on. 

Each of these individual well areas is now divided 
into radial segments, or possibly into radial and 
linear segments. The radius of the segment at the 
injection well, re;, can usually be taken as half the 
distance to the nearest producing well. A value of 
Tep can then be calculated for each well so that the 
sum of the areas of the two radial sectors for each 
well is equal to the area assigned to that well. 

The flood front is then traced through these sec- 
tors in a manner similar to that described for the 
five-spot. In the flood pattern of Fig. 6, for example, 
when the flood front is between the injection well 
and the radial boundary, re; (and prior to break- 
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FIG.6—DIVISION OF IRREGULAR PATTERN INTO 
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through), the injectivity of the pattern referred to 
the difference in pressure between the injection 
well and producing Well 1) is 
or , ii cee 
as oo Twi ~ '¢ 
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Referred to the difference in pressure between the 

injection well und producing Well j, the injectivity 

of the pattern during the same period of the flood is 
| Lh 
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, res B/D-psi 
. (14) 


= On it 6; 4;, reservoir B/D, 
number of producing wells in the pattern ele- 
ment (n = 4 in this example), 
vertex angle of the pattern boundary at pro- 
ducing Well j, radians, 
= total producing rate of Well j,reservoir B/D, 
Est , ft, and 
7 
total area of the isolated pattern element, 
sq ft. 
When the flood front is between rep and the pro- 
ducing well, the injectivity is 


Fis = 
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After breakthrough occurs, the injectivity is 
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where the radius rj; after breakthrough can be as- 
sumed to be constant and equal to the value obtained 
from Eq. 16 using the breakthrough value of areal 
sweep, E,. The angle open to flow from the swept 
region, aj, can be calculated with Eq. 10. 

The isolated element of a producing pattern can 
also be distributed among the producing wells by 
using the differences in pressure between the injec- 
tion well and the producing wells as the criteria 
for assigning the radial-flow sectors to each pro- 
ducing well. Eq. 14 can be written for each produc- 
ing well and the resulting equations solved for the 
appropriate values of rep; according to the differ- 
ences in pressure desired between the injection 
well and the producing wells. One method of doing 
this is by a simple iterative scheme. Eq. 14 can 
be written for each producing well and the results 
combined to give the total injection rate for the 
pattern as 


0.003076 hay ,) @;(P; ~ Pp) 
Qn 








ais SS . ee ee (14a) 

We begin by guessing the ratio of producing rates 
for all neighboring producing wells and calculating 
a set of rep;’s as before. From Eq. 14a an approxi- 
mate total injectivity is calculated. Substituting 


this injectivity into Eq, 14 gives an approximate 
producing rate for each producing well. The ratios 
of these producing rates can be used to re-divide 
the pattern and calculate a better approximation of 
the total injection rate. The initial guess is based 
only on the ratio of producing rates for neighboring 
producing wells and not on absolute rates. When 
the theoretical equations or model data for initial 
injectivities are available for unity mobility ratio, 
the actual injection rate (rather than the rate cal- 
culated from Eq. 14a) is used with Eq. 14 to re- 
divide the pattern area. 





For comparison with experimental model data, 
the injectivity of a nine-spot pattern was calculated 
for an adverse mobility ratio of 10:1, and corner 
wells and side wells producing at equal rates. 
Although the nine-spot pattern generally is not con- 
sidered to be an irregular pattern, it represents a 
system in which the producing wells are not all 
equidistant from the injection well. The injectivity 
was calculated from Eqs. 10, 14, 15, 16 and 17 by 
approximating the pattern area using two different 
criteria for determining the radii of the radial-flow 
sectors—(1) producing rates of the producing wells, 
and (2) the initial injectivities from the theoretical 
equations in Table 1. 

The approximate calculations for the nine-spot 
are shown in Fig. 7, together with porous-model 
results’ for the same conditions. The calculated 
injectivities show very close agreement with the 
model data. The disagreement at early areal sweep 
could be accounted for by wellbore damage in con- 
struction of the model. The slight disagreement in 
the two sets of model data is probably due to in- 
homogeneity in the model and extent of wellbore 
damage at the two corner wells of the model shown 
in Fig. 7. 

BOUNDARY PATTERNS 

To estimate the injectivities of boundary patterns, 
the effect of flow from outside the normal well net- 
work must be taken into account. A study of the 
boundary patterns listed in Table 1 showed that 
the effect of the position of the boundary on the 
injectivity of the outer row of wells approaches a 
limiting value when the boundary is more than one- 
half the well spacing beyond the outer row of wells, 
and it has almost no effect on wells two or more 
rows in from the boundary. All wells except the 
outer row thus can be treated as interior wells, and 
an estimate of injectivity in boundary patterns can 
be made by calculating an effective drainage angle 
6 for the outer row of wells and then tracing the 
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front advance as was done for the regular five-spot 
pattern. For example, if the five-spot pattern with 
all wells producing at the same rate were a boundary 
pattern as shown in Fig. 8, the injectivity of this 
boundary pattern could be estimated as follows. 
The area of this pattern inside the well network 
can be approximated by radial-flow segments iden- 
tical to those of Fig. 2. With these values of ro; 
and rep, the assumed ‘symmetry pattern is then 
represented by the radial sectors shown in Fig. 9, 
Describing the flow in each sector by the radial- 
flow equation and combining gives an initial injec- 
tivity for this assumed-boundary symmetry element 
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res B/D-psi 
ee ee ee . (18) 
producing rate at the boundary well, 

reservoir B/D, 

i = total injection rate of pattern, reservoir 
B/D, and 

6 = effective drainage angle of boundary 
well, radians. 

With all wells producing at the same rate, q) = 
(2/3)i. With the injectivity by the approximate Eq. 
18 set equal to the injectivity by the analytical 
equation in Table 1, the effective drainage angle 6 
can be calculated. For d/r,, = 160, the effective 
drainage angle 6 in Fig. 9 is 2.05 radians. This 
effective drainage angle can be used in the radial 
flow equation, and the injectivity of this boundary 
pattern then can be estimated by tracing the front 
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through these radial sectors in the same manner as 
was done for the interior five-spot pattern. 

Based on the difference in pressure between the 
injection well and producing Well j, the injectivity 
of the five-spot boundary pattern of Fig. 9, prior to 
breakthrough and while the flood front is between 
the injection well and the boundary rg;, is 
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total area within the normal well net- 
work, sq ft, and 

total angle open to radial flow at Well 
j, radians. 
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While the flood front is between the boundary rep; 
and the producing Well j, the injectivity is 
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After breakthrough occurs, the injectivity is 
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where 
fj 8jt 
a + MUI- Hf)... . 3) 
and ry; 1s assumed to remain constant after break- 
een and can be calculated from Eq. 21 using 
the value of E, at the time breakthrough occurs at 
Well j. 

Eqs. 19 to 23 are the general equations for esti- 
mating injectivity in complex well networks and 
reduce to the identical equations of the interior 
flood patterns given in previous sections for the 
five-spot infinite array and for the irregular interior 
patterns. 

A comparison of the approximate injectivity of 
the boundary five-spot pattern of Fig. 9 with the 
injectivity of the interior five-spot pattern of Fig. 2 
is shown in Fig. 10 for a mobility ratio of 10. 


CONCLUSIONS 


When the mobility of the fluid injected into a 
reservoir is different from the mobility of the reser- 
voir fluid, the injectivity changes rapidly during 
the early part of the operation — and then only grad- 
ually until breakthrough occurs. Since most of the 
change in injectivity occurs during the early part of 
the flood while the flood front is nearly radial, reli- 
able estimates of injectivity can be obtained by 
approximating the flood pattern with radial- and 
linear-flow segments. 

Injectivity can be estimated after breakthrough 
by assuming that the flood front can be approximated 
by simple shapes and that the angle the swept region 
makes at the producing well is proportional to the 
rate of flow from the swept region. An estimate of 
breakthrough sweep efficiency is required to calcu- 
late injectivity after breakthrough. 
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The approximate equations of injectivity can be 
applied to irregular interior patterns and to bound- 
ary patterns. 
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NOMENCLATURE 


distance between injection and producing 
wells, ft 
fraction of producing stream which is from the 
swept region 
net pay thickness, ft 
injection rate, reservoir B/D 
specific permeability, md 
pressure, psi 
= producing rate, reservoir B/D 
= radius, ft 
fraction of reservoir swept by injection fluid 
= ratio of mobilities of swept region to unswept 
region 
ratio of producing rates of corner well to side 
well for nine-spot producing pattern 
angle open to flow from the swept region, 
radians 
vertex angle of pattern boundary at producing 
well, radians 
= fluid mobility (k/p), md/cp 
= fluid viscosity, cp 
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= boundary well 
external boundary 
= external boundary of radial sector around in- 
jection well 
external boundary of radial sector around pro- 
ducing well 
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FIG. 1la—~STEADY -STATE PRESSURE DISTRIBUTION 


IN NINE-SPOT PATTERN FOR UNITY MOBIL- 
ITY RATIO AND CORNER-TO-SIDE PRODUC- 
ING RATIO = 0.5. (PER CENT OF TOTAL 
DIFFERENCE IN PRESSURE BETWEEN IN- 
JECTION WELL AND SIDE PRODUCING 
WELL.) d/r,, = 1,000. 
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FIG.11b—STEADY-STATE PRESSURE DISTRIBUTION 


IN NINE-SPOT PATTERN FQR UNITY MOBIL- 
ITY RATIO AND CORNER-TO-SIDE PRODUC- 
ING RATIO = 1. (PER CENT OF TOTAL DIF- 
FERENCE IN PRESSURE BETWEEN INJEC- 
TION WELL AND SIDE PRODUCING WELL.) 
d/r ,, = 1,000. 
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FIG. 1lc-STEADY-STATE PRESSURE DISTRIBUTION IN 
NINE-SPOT PATTERN FOR UNITY MOBILITY 
RATIO AND CORNER-TO-SIDE PRODUCING 
RATIO = 2. (PER CENT OF TOTAL DIFFER- 
ENCE IN PRESSURE BETWEEN INJECTION 
WELL AND SIDE WELL.) d/r,,, = 1,000. 


carried on at the laboratory of The Atlantic Refining 
Co., Dallas. 


8, Muskat, M.: ‘‘The Theory of Nine-Spot Flooding Net- 
works’’, Prod. Monthly (March, 1948) 14. 


APPENDIX 


STEADY-STATE PRESSURE DISTRIBUTION 
IN THE NINE-SPOT PATTERN 


The steady-state pressure distribution in the nine- 
spot pattern for unity mobility ratio was calculated 
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after the method of Muskat!/8 for ratios of corner- 
well to side-well production of 1/2, 1 and 2. For 
the infinite well array shown in Fig. 7, the steady- 
state pressure distribution is given by 


My, cosh 


1(y-kd) ‘ 
d 


/ ™(y-kd) ‘ 


H, =cosh cos —— 


m=+@ 


4rdAhp\x,y) | bly tt 
R y 2m 
/ 2+A 


m=-@® 


/ 
2m +\ 





rye 
(40) 
m=+@ 


= | Resarge pommel 
m=-@® Hom 


. (24) 


RP 


where R =ratio of corner-well to side-well producing 
rates, ‘ 
The calculated pressure distributions are shown 
in Fig. 11 as percentages of the total pressure dif- 
ference between the injection well and the side well. 
The singular point (discontinuity of pressure 
gradient) between the side well and the corner well 
was obtained by cross-plotting of the calculated 
pressure distribution and is only approximate. How- 
ever, the large extent of the near-radial pressure 
distribution is readily apparent. In the case of R = 
1 in Fig. 11, for example, more than 95 per cent of 
the pressure difference between the injection well 
and the corner well could adequately be described 
as radial. week 





Analysis of Gas-Cap or Dissolved-Gas Drive Reservoirs 
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ABSTRACT 


A numerical method of solving the partial differ- 
ential equations which describe the one-dimensional 
displacement of oil by gas has been presented. 
Possible extension of the method to treat multidi- 
mensional flow is discussed, and the limitations of 
this extension are indicated. Using this method, it 
is possible to allow for the existence of a gas cap, 
the presence of any number of gas-injection or oil- 
production wells and the evolution of dissolved gas 
from the oil. It is also possible to allow for variation 
in the cross-sectional area, elevation, porosity and 
permeability of the reservoir. The influence of 
relative permeability and the force of gravity in the 
direction of flow upon the displacement is considered. 

The influence of capillary pressure upon the flow 
and the effect of gravity in the direction perpendicular 
to flow are neglected. The physical properties of the 
fluids are considered to be functions of pressure only, 
and equilibrium between contiguous phases is as- 
sumed. 

The numerical calculations can be readily carried 
out by the use of a digital computer. Several example 
analyses have been performed using the IBM 704 
computer, and about one-third of an hour of comput- 
ing time was required per case. Reservoir behavior 
predicted by use of this numerical method was 
compared to data obtained by other methods for three 
cases — complete pressure maintenance, dissolved- 
gas drive and gas-cap drive. The independent solu- 
tions to these problems were obtained by analytical 
solution, laboratory experiment and field data, re- 
spectively. Agreement of the numerical solution with 
data from these sources was good; this agreement 
establishes the convergence and accuracy of the 
numerical method. 


INTRODUCTION 


Most petroleum reservoirs.can be produced by any 
one of several alternative programs. When a reservoir 
is produced by primary methods, production eco- 
nomics can be influenced by controlling the number 
and location of wells and the flow rate of each well. 
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An even greater influence may be achieved by 
augmenting the recovery of oil obtainable by primary 
methods. This can be accomplished by injection 
of fluids such as water, natural or enriched gas or 
a bank of light liquid hydrocarbons. Selection of 
the most desirable operation requires a means of 
predicting the reservoir behavior which will result 
from each of the several alternative programs. The 
purpose of this paper is to present a mathematical 
method for predicting the behavior of reservoirs 
produced by gas-cap drive, dissolved-gas drive or 
pressure maintenance by gas injection. 

The method described herein takes cognizance 
of phase changes caused by a decline or an increase 
(due to gas injection) in reservoir pressure, of the 
presence of a gas cap and of the effect of gravity 
on the flow of gas and oil. Relative permeability 
relationships are used to define the flow properties 
of the rock. Allowance is made for variation in 
cross-sectional area, elevation, permeability and 
porosity of the reservoir. Both the influence of 
capillary pressure upon the flow and pressure 
gradients in the gas cap are neglected. Whenever a 
liquid phase and a gas phase are in contact, they 
are assumed to be in equilibrium. The physical 
properties of the fluids are considered to be func- 
tions of pressure only. Therefore, if the method is 
to be used to predict the effects of a gas-injection 
program, mixtures of the injected gas and formation 
crude should-have the same physical properties as 
mixtures of formation gas and crude. 

The equations to be presented in this paper apply 
only to a one-dimensional case; therefore, they 
neglect the influence of gravity in the direction 
transverse to the flow. As is well known, this 
gravitational influence may lead to overriding of 
oil by gas. Consequently, this procedure as presentec 
is most applicable to long, thin reservoirs for which 
gravity overriding is not important. On the other 
hand, the equations presented can be generalized 
to treat multidimensional flow and, hence, to consider 
gravity overriding, if desired. A word of caution 
on two points is advisable here, however. First, 
the authors have not demonstrated the accuracy of 
the numerical technique for multidimensional flow. 
Second, and more important, capillary pressure will 
often be of importance in multidimensional problems. 
Obviously, in such cases a generalization of the 
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present method would be inadequate. 
These limitations of the present method are less 
restrictive than limitations on previously published 
methods of predicting gas-drive reservoir behavior. !~9 
For example, the method of Ref. 1 is restricted by 
the assumption that pressures and saturations are 
uniform throughout the reservoir, and the method of 
Ref. 2 is limited to dissolved-gas-drive reservoirs 
involving a single well. Ref. 3 through 8 are suitable 
only for complete pressure maintenance by gas 
injection. Ref. 9 presents a method capable of 
handling pressure decline, but the mathematical 
relations used are not rigorous, the method of 
computation is tedious and it does not provide a 
way of calculating pressure distribution throughout 
the reservoir. 

The method presented in the paper differs from 
that of Ref. 9 by being based on solution of a more 
rigorous mathematical statement of the problem. 
This statement consists of two simultaneous partial 
differential equations and suitable boundary condi- 
tions. These equations are solved by a systematic 
approximation technique which has been adapted to 
solution by electronic computation. This relieves 
the engineer of most of the labor involved and makes 
feasible the completion of a reservoir analysis in 
approximately 20 minutes on an IBM 704 computer. 
The calculations provide complete pressure and 
saturation profiles at frequent time intervals, as 
well as gas-oil ratios and cumulative production 
data for each well. 

Derivation of the differential equations stating 
the problem is presented first. This is followed 
by transformation of the equations into finite- 
difference equations. Finally, the convergence and 
accuracy of the method is evaluated by comparing 
computed results with results obtained from analyt- 
ical solution, laboratory experiment and field 
observation. 


THE DIFFERENTIAL SYSTEM 


The fundamental equations used to define the 
reservoir behavior are the equations of continuity 
for gas and stock-tank oil and Darcy’s law of flow 
for each phase. The continuity equation for gas 
allows for the presence of a distributed sink for 
fluid production of strength V, for oil and V, for 
gas. This equation is 


0 
¢ aye + u,AB,R,) -V,R, - V, = 


pA <(s,5 * OBR) 6.6.9.0 60040 6nis (1) 


where 5 =the standard cubic feet of gas per cubic 

feet of gas at reservoir conditions, 

Vp =the rate of free-gas production per unit 
length of the reservoir, and 

Vo =the rate of oil production per unit length 


_ of the reservoir. ; ; 
Eq. 2 is the continuity equation for oil. In this 





lReferences given at end of paper. 
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expression, the volatility of the stock-tank oil has 
been neglected. 


- © (4 AB) -V, = $A 2 (s.B,) -oe @ 
Ox at 


Darcy’s law of flow for each phase, 





kk, [aP a 
u,=- - tI Lae ave (3) 
He \de Oz 
and 
kk, { aP a 
u,=-—* —+e.9 — esrees (4) 
He Oz Oz 


where z = the elevation of the reservoir above an 
arbitrary datum plane, and the equation, 


complete the statement of the differential problem 
except for boundary and initial conditions. 


Eqs. 3, 4 and 5 may be used to eliminate uo, ug 
and S, from Eqs. 1 and 2. This would result in two 
differential equations containing two dependent 
variables, S$, and P. These two equations may be 
solved simultaneously to obtain the desired solution. 
Instead of taking this approach, it is desirable first 
to combine Eqs. 1 and 2 to eliminate time derivatives 
of saturation and, then, to use the new equation and 
Eqs. 2, 3, 4 and 5 to define the problem. The advantage 
of doing this is that the new equation can be approx- 
imated by a difference equation having only pressures 
at an advanced time level as unknowns. This gives 
rise to a set of equations which can be solved to 
obtain a complete pressure distribution at the new 
time. These new pressures canbe used ina difference 
approximation of Eq. 2 to yield saturations at the 
new time. Thus, the solution of the two equations 
is divided into two steps, solution of a pressure 
equation and solution of a saturation equation. The 
necessity of solving for pressures and saturations 
simultaneously is avoided by this procedure. 


The desired combination of Eqs. 1 and 2 is effected 
by multiplying Eq. 1 by B,, multiplying Eq. 2 by 
(5 — BoRs), replacing Sg by (1 — Swe — Sg), expand- 
ing the product derivatives on the right-hand side of 
the equation and replacing terms of the form 





aC(P) , aC(P) OP _ 6, dP 
at dP at at 


Here C or C(P) indicates a function of pressure 


only, and C” denotes a, Adding the resulting 
dP 


equations gives 















a 
- By = (UgAB + UAB Re) ~ (5 - BR.) 


a 
= (UeAB,) - B,(V,R, + V,) - (8 - BL RV, 
x 


= pAIS,[B, (5 - B, R,)’— Bi(5 - B, R,)| + 


[1 - S,,,] [B,(B, R,)’+ By(8 - B, R, I} e - 
t 
PALD(P, S,)] ~ Pe are Aaa akan on (6) 
t 


where the function, D(P, S,) is defined by the last 
equality. 

Next, uo, ug and So are eliminated from Eqs. 6 
and 2 by use of Eqs. 3, 4 and 5, which results in 
the following equations. 


3 kk, »Ad oP De 
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All that is needed in addition to Eqs. 7 and 8 to 
define the reservoir behavior is a statement of 
suitable boundary and initial conditions. 

The initial conditions in most gas-cap-drive 
situations are that the gas saturation throughout 
the reservoir is zero and that the pressure at the 
original gas-oil contact is specified. Since the 
reservoir fluids are initially at rest, the pressure 
gradient at any point in the reservoir is equal to the 
density of the oil times the elevation gradient at 
that point. This completely defines the initial pres- 
sure distribution in the reservoir. Stated symboli- 
cally, these initial conditions are 


S 


0 for O<2<X, t=0; 





& 
P=P,forrz2=0, ¢=0 ; 
and a 
dP(@) 1 splig ae Me 0<2<X,t=0 
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where X = total length of the reservoir, and 
P, = initial gas-cap pressure. 

The boundary condition at the original gas-oil 
contact is that the rate of gas entry at that point, 
ig, is determined by the expansion of the gas-cap 
gas which has a volume V. The amount of gas in the 
gas Cap at any pressure is V § (P), and the rate at 
which gas enters the oil zone is —V 5° OP 

Ot 

At any production-well location, it is desirable to 
specify the rate of stock-tank oil production qg 
provided this rate does not exceed the productivity 
of the well. This oil rate is related to the gas- 
production rate through the relative permeabilities 
and fluid properties existing at that time at the 
location in question. The quantitative relationship 
is derived from material balances on the differential 
element Ax of the reservoir which contains a pro- 
duction well. A material balance on the oil phase 


yields 
B,Alu,), = 9, +B, Alu), 


= Lg 


where subscript 1 indicates conditions at x — \x/2 
and subscript 2 conditions at x + Ax/2. 
The gas-phase material balance yields 


bA(u,) | = 9,+ 5A(u,), 
Darcy’s law of flow is used to eliminate (ug),, 


(ug), (uo), and (uo), from these equations; then the 
resulting equations are combined to eliminate 


oP _ si 
i, Wael, 
seein (= | “(): 


6 keg Fo 
Oo B, k 


the result is 


hy +d, cee Scennnnncannes (9) 


ro Hg 


If the well is located at the end of a reservoir, 
then (uo)2 = (ug)2 = 0. In this case, Eq. 10 is the 


result of this analysis and should be used in place 
of Eq. 9. 





q Ok, op kk, ,A8g a 
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Note that in both Eqs. 9 and 10 q, represents only 
gas produced as free gas. A quantity of dissolved 
gas equal to q, R, is also produced with the oil. 


If the productivity of the well is less than the 
desired production rate, the boundary condition 
becomes one of constant pressure at the well. 
Since the one-dimensional analysis presented in 
this paper neglects the pressure drop due to radial 
flow around the well, the constant pressure employed 
in the calculations should be greater than the 
actual bottom-hole pressure to be maintained in the 
well. The difference in these two pressures should 
equal the pressure drop caused by radial flow into 
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the well and can be estimated from radial flow 
formulas. 

At the reservoir boundaries, there must be no 
flow of either gas or oil. This is accomplished by 
setting Uo = Ug = 0. 

At injection wells (if there are any), the rate of 
gas injection may be arbitrarily specified. The gas 
saturation at gas-injection wells and at the gas- 
cap contact is set equal to that saturation existing 
when the rock contains residual oil. 

Specifying the boundary conditions completes the 
statement of the problem in differential form. The 
next step is to write suitable difference approxima- 
tions to Eqs. 7 and 8, 


THE DIF FERENCE SYSTEM 


There exist many difference equations which under 
casual investigation appear to be valid analogues 
of Eqs. 7 and 8. Some of these provide adequate 
approximation, while others do not. Since there is no 
rigorous analysis to demonstrate that the solution 
of these analogues converges to that of the differ- 
ential system, the final test must be experimental 
computing. Analysis of similar, but simpler, systems 
serves as a useful guide in selecting equations to 
be tested. 

The set of difference approximations to Eqs. 7 
and 8 presented in this section have been subjected 
to extensive testing, and the results are discussed 
in the following section. This testing has shown 
the set to be suitable for the prediction of reservoir 
behavior for gas-cap drive, dissolved-gas drive and 
pressure maintenance by gas injection. However, no 
claim is made that this is the only set of equations 
which will satisfactorily approximate Eqs. 7 and 8. 

In writing finite-difference approximations to 
Eqs. 7 and 8, the distance axis is divided into J 
increments Ax in length, and the time scale is 
divided into increments At in length. The following 
definitions are used. 

x; =] Ax, 

X = J Ax = total length of reservoir, 

t, =n \t, and 

Pin ~ P(x;, t, ), etc. 

Derivatives of the form, 


> 
2 lovpye(s,) = 
Ox Ox 


to be evaluated at x; and t,,, are replaced by terms 
of the form, 


Oxy (Fat, ntt ny) ale FF net F-1, nd) 
a 2 ot“ j, A ia a, *j-l, ‘ 
j+%,n™j,n j-%,njel,n ) 





Az 


in arriving at difference approximations of Eqs. 7 
and 8. 

Attention is called to the fact that the coefficient 
E, which is a function of saturation only, is evaluated 
at the upstream end of the interval. That is, if the 





JUNE, 1961 


per cent 


Relative Permeability, 


Gas Saturation, fractional 


pressure difference involves points j and j + 1, E is 
evaluated at point j. This is the upstream end of the 
interval since in the problem being considered flow 
is from point 7 to 7 + 1. If E were evaluated at the 
midpoint of the interval, like C (P), then the 
numerical procedure would not give solutions which 
converge to a true solution of the differential sys- 
tem. If the equations presented herein are generalized 
to treat multidimensional flow, care must be taken to 
evaluate the function E at the upstream end of the 
interval whenever it appears in approximations of 
distance derivatives. Otherwise, the extension of 
100 
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FIG. 1 — RELATIVE PERMEABILITY-SATURATION RE- 
LATION, PRESSURE MAINTENANCE. 
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FIG. 2 — SATURATION PROFILE, PRESSURE MAIN- 
TENANCE. 









the procedure to two or three dimensions is perfectly 
straight-forward. 


Time derivatives such as Go are approximated 
j 
P, 


as jp, n+i1 


At 
mated as distributed sinks of strength (2) for 
Ax jm 





i} ® Production wells are approxi- 


. q ae 
oil and —? for gas. Injection wells are treated as 
[Nx . 


: ae. 
Ax 
source and sink terms satisfies the boundary con- 
ditions of specified rates at the wells. 

In writing the difference equations, subscripts 
were dropped from Bo, Rs and Sg for simplicity. 
The approximation to Eq. 7 is implicit in pressures 
and is given by 


sources of strength . The inclusion of these 
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FIG. 3 — RELATIVE PERMEABILITY-SATURATION RE- 
LATION, DISSOLVED-GAS DRIVE. 
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FIG.4 — EXPERIMENTAL PRODUCTION HISTORY, 
USED AS BOUNDARY CONDITION TO COMPUTE 
BEHAVIOR OF DISSOLVED-GAS DRIVE. 
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where ig, Yo, 4g > 0. 

There are (J + 1) grid points in the reservoir 
corresponding to 7 =0, 1,2... J-1, J. The reser- 
voir may be divided into segments of length \x, 
centered about grid points. Eq. 11 as applied to the 
point 7 represents a combination of the two material- 
balance equations for the segment centered on }. 
At the end points j = 0 and j = J, there exists only 
one-half of a segment; hence, the right-hand side of 
Eq. 11 should be divided by two when applying it 
to these points. 
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Maintaining the proper boundary conditions at the 
points j = 0 and j = J also requires that k — 1 and 
k, be set equal to zero in Eq. 11. The resulting 
form of the equation places a restriction on the 
pressure such that the rate of flow at these points 
is equal to that specified by the i, qo, and qx 
terms of Eq. 11. If no production or injection occurs 
at j=], then ig =q,=q,=0 and the no-flow boundary 
conditions of up =ug=0 are satisfied. The boundary 
condition at the gas-cap contact is satisfied by 
setting ig at j = 0 equal to 


(P 


o,n+1 F ow 
At 
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As stated earlier examination of Eq. 11 reveals 


that, if conditions of saturation and pressure are. 


known at time level t,, the only unknowns in Eq. 11 
are pressures at time ¢,,,;. The equation may, there- 
fore, be put in the form, 


aP iy nti + Pi ng + CP itt ned = 9; 
where a;, b;, cj and d; are known constants for each 
point 7. The equation for each point involves three 
unknowns, except for the first(j7=0) and last (j= J) 
equations which contain only two. There are, al- 
together, J] equations and J unknowns. This system 
of linear simultaneous equations can readily be 
solved by a method described in the literature!® to 
yield a set of pressures at the new time level ty44. 
These new pressures are used in a difference 
approximation of Eq. 8 to find the new saturation at 
tn41- The approximation used for Eq. 8 is 
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At 7 = 0, the original gas-oil contact, the gas 
saturation is held constant at a value corresponding 
to the residual oil saturation. At j = J, as for Eq. 
11, it must be noted that k; = 0 and that only one- 
half of a segment is present. 

When the pressure at a well is specified rather 






JUNE, 1961 








than the flow rate, this pressure is first corrected 
to allow for the pressure drop due to radial flow 
into the well. Then this corrected pressure is used 
in lieu of Eq. 11 at the grid point corresponding to 
the well. Writing Eq. 11 for all points except those 
at which pressures are specified and combining 
this set of equations with the specified pressures 
results in a system of equations which can be 
solved to yield a complete pressure distribution at 
time level n + 1 (assuming pressures and saturations 
were known at time level n). Then Eq. 11 written 
for the points at which the pressure was arbitrarily 
specified will contain two unknowns, qo and qg, 
and can be combined with either Eq. 9 or Eq. 10, 
as appropriate, to yield qo and qg for the well. 
Eq. 12 is solved explicitly for the saturation at 
time level n + 1. As with all explicit equations, 
there is an upper limit on the size of the time step 
that can be employed. This limit cannot be estab- 
lished precisely for. the general case, but for the 
case of constant coefficients (i.e., the pressure- 
maintenance problem where fluid properties are 
considered invariant), this limit is given by 


Ard 


At< |———____ 
(u, + uf 


minimum 


The maximum value of (uw, + ug) is employed in 

dig 
dS 
evaluated at the gas saturation existing at the gas 
front. 

This limitation was derived from a heuristic 
stability analysis and verified by experimental 
computing for the case of pressure maintenance. 
While it was found that a slightly less stringent 
time-step limitation could be applied for gas-cap or 
dissolved-gas-drive calculations, this equation 
serves as a useful guide in selecting the time step. 


this equation. The quantity le is equal to 


MATERIAL BALANCE EQUATIONS 


The difference equations presented in the preced- 
ing section are consistent with the conservation of 
matter. Thus, the amount of gas and oil should be 
conserved during the course of the calculations. It 
is advantageous to monitor the accuracy of the cal- 
culations by means of an over-all material balance 
on both gas and oil, since such balances represent 
independent checks on the accuracy with which the 
difference equations are being solved. In this way, 
machine or other errors may be detected. 

To compute an over-all material balance it is 
necessary to integrate to obtain both the quantity of 
oil and the quantity of gas present in the reservoir 
as a function of time. It also is necessary to inte- 
grate to obtain the cumulative production of gas and 
oil. The integrals are of an elementary form and 
will not be repeated here. However, it should be 
noted that, to make the material balance exact 
except for machine round-off error, these integrals 
must be evaluated by the trapezoidal rule. 
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EVALUATION OF THE METHOD 


The method described in the preceding sections 
is intended for application to a wide variety of gas- 
drive problems including dissolved-gas drive, gas- 
cap drive, complete pressure maintenance and 
pressure-buildup operations. Therefore, in evaluating 
the validity of the method, it was desirable to solve 
a number of problems by the numerical method of 
calculation and to compare the answers to ones 
obtained by analytical methods, laboratory experi- 
ment or field data. Three problems were selected 
— complete pressure maintenance, dissolved-gas 
drive and gas-cap drive. In all three cases, gravity 
was an important factor in the recovery mechanism. 
The first problem, complete pressure maintenance, 
was solved analytically by the Buckley-Leverett 
technique. Laboratory data were obtained for a 
dissolved-gas-drive reservoir, and field data were 
available for a gas-cap-drive reservoir. 

The numerical calculations were programed for 
the IBM 704 in machine language using the SHARE 
assembly program. The total machine storage re- 
quirement for input-output routines, program and 
data storage is approximately 7,000 words. If the 
reservoir is divided into 20 grid intervals, which 
is adequate for the problems described, the compu- 
tations require six seconds per time step. About 
200 time steps are required in a typical case, and 
this results in a total computing time of approxi- 
mately 20 minutes per case. Increasing the number 
of grid intervals increases the number of time steps 


required, so that total computing time is proportional 
to the square of the number of intervals used. 


COMPLETE PRESSURE MAINTENANCE 


The results of complete pressure maintenance 
were computed for a hypothetical reservoir 10,000-ft 
long, 25,000 sq ft in cross-sectional area and with 
an angle of inclination of 6° from the horizontal 
direction. The reservoir was initially 75 per cent 
saturated with oil; the other 25 per cent saturation 
was connate water. Production was taken only from 
the bottom extremity of the reservoir, and gas was 
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Pressure at Top of Column, psig 
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FIG. 5 — PRESSURE BEHAVIOR AT TOP OF COLUMN, 
DISSOLVED-GAS DRIVE. 





injected into the upper extremity at a rate of 500,000 
scf/day. This rate of injection just balanced the 
total rate of production so that the pressure at 
the lower end of the reservoir was constant at 
1,000 psig. Reservoir permeability was 1 darcy, and 
the porosity was 25 per cent. The relative permea- 
bility-saturation relationship for the reservoir matrix 
is shown in Fig. 1. 

This problem was solved under the assumption 
that the pressure drop in the reservoir was small 
enough so that the fluid properties could be treated 
as invariant. The fluid properties employed in the 
calculations were: oil viscosity, 0.5 cp; gas vis- 
cosity, 0.014 cp; oil density, 45.5 lb/cu ft; gas 
density, 0.5 lb/cu ft; oil formation volume factor, 
0.8; gas solubility, 600 scf/STB; and gas expansion 
factor, 82.5 scf/reservoir cu ft. 

The saturation profile after 8.2 years of injection, 
as calculated by the Buckley-Leverett4:5 analytical 
technique, is shown as the solid curve in Fig. 2. 
Results of the numerical calculations are also shown 
in this figure. The open circles were computed 
with Ax = 250 ft, the squares with Ax = 500 ft and 
the solid circles with Ax = 1,000 ft. The ratio At/ 
Ax was equal to 0.41 x 10-3 years/ft in all three 
calculations. Behind the gas front, the data from 
all three numerical calculations give good agree- 
ment with the analytical solution, but at the front 
it is only fair for Ax equal to 500 or 1,000 ft. With 
Ax equal to 250 ft, the computed front is relatively 
accurate; further refinement of the size of Ax, 
while keeping At/Ax constant, would further in- 
crease the accuracy. 

Refinement of At while keeping Ax constant at 
500 ft did not appreciably alter the solution shown 
in Fig. 2. 

Application of the numerical technique to this 
problem constitutes a severe test of the method 
because of the saturation discontinuity which devel- 
ops during the displacement. Nevertheless, the 
data given in Fig. 2 show that the approximate 
solution can be made to approach the analytical 
solution as closely as desired by refining the size 
of Ax. 

4 
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Oil Phase Potential Difference, psig 
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FIG. 6 — VARIATION OF THE DIFFERENCE IN OIL- 
PHASE POTENTIAL AT TOP AND BOTTOM 
OF COLUMN, DISSOLVED-GAS DRIVE. 
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TABLE 1 — PHYSICAL PROPERTIES OF FLUIDS FOR 
DISSOLVED-GAS-DRIVE SYSTEM, TEMPERATURE = 75°F 
fa Oil 

hie Solu- Formation 
sure bility* Volume 
(psig) Factor 

0 0 1.000 1.00 1.55 0.115 

5 6.26 0.991 - - a 

10 11.73 0.980 1.70 1.31 0.198 

15 17.48 0.967 - a ee 

20 23.58 0.950 2.43 1.11 0.280 

25 30.32 0.930 - ss as 

30 38.00 0.908 3.17 0.93 0.364 

35 47.62 0.882 ~ ‘a ~ 

40 62.06 0.854 3.94 0.76 0.450 


*Cubic feet of gas at stondard conditions per cubic foot of 
stock-tank oil. 


Expansi- Oil Gas 
bility Viscosity Density 
of Gas** (cp) (Ib/cu ft) 








** Cubic feet at standard conditions per cubic foot at reservoir 
conditions. 





DISSOLVED-GAS DRIVE 


Dissolved-gas-drive calculations were used to 
simulate the behavior of a laboratory model. This 
model consisted of a rectangular column packed 
with sand and mounted with the long dimension 
vertical. The dimensions of the model were 3/8 x 
1/2 x 72 in. It was packed with 1/2 Wausau sand 
to a porosity of 44.1] per cent and a total permeability 
of 13.0 darcies. Relative permeability data for the 
sand are presented in Fig. 3. These data were 
determined experimentally by the displacement 
technique of Johnson, et al.41 

The oil phase charged to the column was a light 
refined hydrocarbon (Bayol D) saturated with propane 
at a temperature of 75°F and a pressure of 36.64 
psig. The experiment was conducted at a constant 
temperature of 75°F; some of the properties of the 
equilibrium phases of Bayol D-propane mixtures at 
this temperature are given in Table 1 as a function 
of pressure. These properties are the oil formation 
volume factor and viscosity, and the gas expan- 
sibility, solubility and density. The gas viscosity 
was constant at 0.008 cp. The oil density varied 
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—— Computed Data 
© Experimental Data 


Free Gas Production, % of original 
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FIG. 7 — FREE-GAS PRODUCTION HISTORY, DIS- 
SOLVED-GAS DRIVE. 


JUNE, 1961 


linearly from 47.5 lb/cu ft at a pressure of 0 psig 
to 45.0 lb/cu ft at a pressure of 40 psig. At the 
start of the experiment, the model was saturated 
with oil containing dissolved gas, there being no 
connate water present. Production of oil was initiated 
from the bottom of the column by reducing the back- 
pressure maintained at this point. The resulting 
cumulative production of stock-tank oil is shown in 
Fig. 4 as a function of elapsed time. 

A material balance applied to the experimental 
data indicated that the oil and gas in the column 
after the initiation of production were not in complete 
equilibrium. This lack of equilibrium probably was 
caused by supersaturation of the liquid phase. The 
average degree of supersaturation during the early 
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FIG. 8 — RELATIVE PERMEABILITY-SATURATION RE- 
LATION, GAS-CAP DRIVE. 








foe] 
oO 


Cross Sectional Area, thous. of sq. ft. 
S 








| | | 
2 4 6 


Distance from Gas Cap, thous. of ft. 


FIG. 9 — VARIATION OF CROSS-SECTIONAL AREA, 
GAS-CAP RESERVOR. 








TABLE 2—PERTINENT RESERVOIR DATA —GAS-CAP DRIVE 


Original and Saturation Pressure (psig). . . » . «+ 2,700 
Reservoir Temperature (°F)... 1s ee we we eee 181 
Volume of Oil Zone (acre-ft), . . . 2 ee eo 20,986 
Volume of Effective Gas Cap (acre-ft) ... . 10,883 
Formation Dip (ft/mile). « « « «2 0 6 0 e ° 226 
Average Permeability (md) . .. +. +s eos : 165 
Average Porosity (per cent). . 2... ee eee ee 23 
Average Connate-Water Content of Oil Zone (per cent) . 28 
Average Connate-Water Content of Gas Zone (per cent). 25 
Original Oil in Place (STB)... . 19,840,000 
Original Gas-Cap Gas in Place (Mscf) - « 13,969,000 





part of the experiment was estimated to be slightly 
in excess of 1 psi. Inconsistencies in the material 
balance on gas of about 4 per cent prevented esti- 
mation of the degree of supersaturation late in the 
experiment. 

In computing the behavior of this model, the rate 
of production of stock-tank oil at any time was 
derived from the data of Fig. 4. This rate schedule 
was entered into the calculations as the boundary 
condition at the production end of the model, and 
the resulting gas production and pressure behavior 
of the model was calculated by the previously de- 
scribed numerical technique. 

Both the results of calculation and experiment 
are shown in Figs. 5, 6 and 7. In all three figures, 
the abscissa is the cumulative quantity of oil 
produced. The ordinate in Fig. 5 is the pressure at 
the top of the column; in Fig. 6 it is the difference 
in the oil-phase potential (P + pg g A z) measured 
at the top and bottom of the column, and in Fig. 7 
it is the cumulative quantity of free (undissolved) 
gas produced. In these figures, the solid curves 
represent the calculated data; the open: circles 
are data points which were experimentally measured. 

The computed pressure behavior at the top of the 
column (Fig. 5) agrees quite well with that observed 
experimentally, especially early in the experiment. 
Agreement was not as good late in the experiment 
when the pressure was declining very rapidly, making 
it experimentally difficult to measure accuratély 
the cumulative quantity of oil produced. Note, 
however, that the final cumulative oil recovery 
measured after the final blowdown to zero pressure 
agreed very well with the calculated ultimate 
recovery. 





Total Oil Production Rate, st. bbi/day 








a 
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Time, yrs. 


FIG. 10 — VARIATION OF TOTAL PRODUCTION RATE, 
GAS-CAP RESERVOIR. 

The difference in oil-phase potential measured 
between the top and bottom of the column (Fig. 6) 
was relatively small, 2.5 to 3.5 psig, and the 
pressure difference was of course smaller, 0.6 to 
1.6 psig. The experimental data exhibit appreciable 
scatter; but, even so, the average deviation of the 
experimental from the computed values is + 12 per 
cent. 

The data of Fig. 7 indicate that the first free-gas 
production occurred when about 15 per cent of the 
original oil had been produced. Thereafter, the 
agreement of the experimental and computed data 
is reasonably good. In fact, the comparisons shown 
on Figs. 5, 6 and 7 indicate that the computed and 
experimental reservoir behavior agree within the 
probable error of the experiment. 


GAS-CAP DRIVE 


Kirby, et al,12 published a limited amount of 
production history for a gas-cap-drive reservoir in 
which there was no active water drive. Data for 
this reservoir are given in Table 2. Table 3 gives 
data on fluid properties as a function of pressure 
for the reservoir fluids. Relative permeability 
data for the reservoir are shown in Fig. 8. These 





TABLE 3 — PHYSICAL PROPERTIES OF FLUIDS FOR GAS-CAP-DRIVE SYSTEM, TEMPERATURE = 181°F , 


Pressure 


Oil Vol. Expansibility 
(psig) ** 


Solubility* Factor of Gas 








0 0 1.000 0.8 


0.930 17.0 
0.888 34.0 
0.865 51.5 
0.845 69.3 
0.823 87.2 
0.802 106.5 
0.782 126.0 
0.760 144.0 
0.738 159.0 


300 17.0 
600 31.3 
900 42.8 
1200 51.0 
1500 63.0 
1800 75.5 
2100 88.5 
2400 102.8 
2700 119.0 


Oil Gas Oil enact 
Viecory ee itveu)—(ltveu td 
0.816 0.0117 50.3 ; 0.03 
0.664 0.0131 47.6 0.75 
0.576 0.0145 46.5 1.48 
0.515 0.0157 45.8 2.20 
0.470 0.0168 45.0 2.95 
0.435 0.0177 44.2. 3.73 
0.407 0.0185 43.5 4.52 
0.383 0.0194 42.7 5.22 
0.360 0.0201 42.0 5.93 
0.338 0.0209 41.2 6.65 


* Cubic feet of gas at standard conditions per cubic foot of stock-tank oil. 
** Cubic feet at standard conditions per cubic foot at reservoir conditions. 
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FIG. 11 — GAS-CAP PRESSURE, GAS-CAP DRIVE. 


data were measured on cores taken from the field. 
Fig. 9 shows the variation of cross-sectional area 
with location in the reservoir, and Fig. 10 shows the 
oil-production schedule for the reservoir. 

All of these data were employed in the numerical 
computation of the reservoir behavior. Producing 
wells were distributed throughout the reservoir as 
indicated by the field map given in the reference 
paper. 

In Fig. 11 computed gas-cap pressures are compared 
to those observed in the field. The solid curve is 
computed, and the circles represent the field data. 
Agreement is quite good; the best curve through the 
experimental points would fall a little above the 
computed curve, with a maximum deviation of 15 psi. 

For the production rates used in this reservoir, 
oil-zone pressures differ only slightly from the 
gas-cap pressure. This is illustrated in Fig. 12, 
which shows shut-in bottom-hole pressures for a 
well 1,900 ft from the gas cap as a function of 
total cumulative production. As in Fig. 11, the 
solid curve represents the computed data, and the 
circles are the observed data. The data presented 
in Fig. 12 were not presented per se by Kirby, 
et al, but were taken from their original worksheets. 
Kirby, et al, chose to present an average pressure 
of all the wells in the field, rather than the pressures 
for a single well. Also, they adjusted the observed 
pressures at all wells to a common elevation by 
subtracting the quantity (po g A z)from the observed 
pressures. Here A z is the difference in elevation 
between the well being observed and the reference 
elevation. The computed pressures in this figure 
are a little higher (15 psi or less) than the computed 
gas-cap pressure at corresponding amounts of 
production. As for the gas-cap pressures, the com- 
puted pressure data are inreasonably good agreement 
with the observed data. 


CONCLUSIONS 
A numerical method of solving the partial differ- 
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FIG. 12 — OIL-ZONE PRESSURE 1,900-FT FROM GAS- 
CAP CONTACT. 


ential equations which describe the displacement 
of oil by gas has been presented. This method 
and the assumptions made in deriving it are described 
in the paper. 

Results presented in the paper lead to the follow- 
ing conclusions. 

1. Comparison of data computed by this method 
to the analytical solution for complete pressure 
maintenance demonstrates the convergence and 
accuracy of the numerical method. 

2. Comparison of computed results to data obtained 
from a laboratory model produced by dissolved-gas 
drive shows the applicability of the method to this 
mode of production. 

3. Finally, a comparison of predicted data with 
field data from a gas-cap-drive reservoir demonstrates 
its applicability to a field problem. 
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DISCUSSION 


J. C. MARTIN 
MEMBER AIME 


The authors have considered a very difficult 
problem and one in which many questions remain to 
be answered. The partial differential equations are 
very complex, and the finite-difference method of 
solution has not been established rigorously. It is 
known that finite-difference methods can fail to 
approximate the true solutions of partial differential 
equations. Thus, there is reason to doubt the validity 
of the method from a mathematical viewpoint. 

The results of reservoir analysis may lead to 
large capital investments and may be a prime factor 
in determining the amount of oil recovery. Thus, 
there is an economic incentive for the results of 
reservoir analysis to reflect the true nature of the 
reservoir. Misleading results can be worse than no 
results at all if they lead to poor‘investments and 
the waste of important national resources. Thus, 
there is a considerable justification for being very 
critical of the validity of new methods of reservoir 
analysis. 

The method of analysis presented by the authors 
involves the solution by finite-difference approxi- 
mations of two simultaneous partial differential 
equations containing the fluid pressure and one of 
the fluid saturations as dependent variables, and 
distance and time as the independent variables. 
Since two dependent variables can vary with dis- 
tance and time, it seems reasonable to expect that 
the authors would present examples where the pres- 
sure and the saturations have considerable varia- 
tions with both time and distance. Also, one would 
expect that numerical results would be shown to 
approach the results obtained by analytical methods 
as the increments in time and distance are reduced. 

In the example of complete pressure maintenance, 
the numerical results do not appear to converge to 
the analytical solution of the partial differential 
equations. The analytical solution is contained in 
the well known paper by Buckley and Leverett. 1 It 
is maintained by some that the Buckley-Leverett 
solution is not the true analytical solution and that 





lReferences given at end of Discussion. 
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the multivalues given by Buckley and Leverett occur 
because they made a mathematical transformation 
which is not valid in the vicinity of the saturation 
discontinuity. Thus, there is some disagreement 
concerning the Buckley-Leverett solution. For this 
reason, it is discussed in the following paragraphs. 

The Buckley-Leverett solution can be obtained 
by well established mathematical methods of the 
type employed in Ref. 2. The partial differential 


equation for the gas saturation is 


ibid 
rs Fane ) — 28(Po ~ Pe) cs J - 3s 
Ox 8 Ho & Ot 


(1) 


where a = the dip angle, 
V(t) = a known function of time, 


F(S,) = — 
1S 


Ho Rp 
G(Sg) = ko/k F(S,). 
Eq. 1 is a first-order partia| differential equation 
of a type for which there is a well known method of 


obtaining solutions (see, for ¢xample, Chapter 4 of 
Ref. 3). The analytical solution of Eq. 1 is 


, and 


x = FS.) a ase 
6 = d Uo 
5 Se Gee ee ab ae Qe 
where W (Sg) is the initial gas saturation. For the 
case considered by the authors where the initial 
gas saturation is zero, W(S,g) is zero. That Eq. 2 
satisfies Eq. 1 can be verified by direct substitu- 
tion. Eq. 2 yields the same multivalue saturations 
profiles as obtained by Buckley-Leverett. Thus, 
the Buckley-Leverett solution is the true analytical 
solution, and it contains multivalued saturations. 
Furthermore, the analytical solution does not con- 
tain a discontinuity. 
In formulating the partial differential equation, it 
is assumed that the saturations are continuous 
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functions of distance. The presence of a discontin- 
uity violates this assumption and across the dis- 
continuity the partial differential equation breaks 
down. If it is assumed that a discontinuity exists 
at a given point and a given time, continuity con- 
siderations yield the relation for the velocity of 
the discontinuity. Buckley-Leverett, followed by a 
number of investigators, have used the procedure 
of first determining the analytical solution and then 
eliminating the presence of multivalue saturations 
by employing discontinuities which satisfy condi- 
tions. 


The results presented in Fig. 2 do not approxi- 
mate the true analytical solution but appear to ap- 
proximate the combination of the analytical solution 
and the discontinuity which satisfies material 
balance. The authors should explain why thé re- 
sults do not approach the true analytical solution 
since the finite-difference relations are represented 
as being approximations to the partial differential 
equations. Furthermore, in the vicinity of the dis- 
continuity the numerical results are a very poor 
approximation to the combination of the analytical 
solution and the discontinuity. Thus, it is not clear 
that the finite-difference method would yield results 


Many of Martin’s comments concern the example 
case of complete pressure maintenance. The ana- 
lytical statement of this problem was formulated by 
Buckley and Leverett! as a first-order quasilinear 


hyperbolic differential equation. This class of 
differential equation often arises in simplified 
theories which ignore mechanisms of dissipation 
such as viscous stresses, heat conduction, ohmic 
loss or, as in the present case, capillary pressure. 
It is the rule, rather than the exception, that the 
physically meaningful solutions of these problems 
develop singularities (discontinuities) after a finite 
time.2 Lax? notes that, while the classical solution 
of the differential equation cannot be continued 
after this time, it is possible to continue it in a 
generalized sense. This kind of generalization is 
dictated by the integral version of the conservation 
law. Lax refers to this generalized solution as a 
weak solution and gives a formal definition of it. It 
develops that these weak solutions are not deter- 
mined uniquely by their initial values. Therefore, 
an additional principle is needed for developing a 
relevant subclass. In this aspect, the problem is 
similar to the solution of a quadratic equation. A 
unique answer cannot be obtained from the equation 
alone; additional information is necessary to reject 
the extraneous root. 

With this background,-we are prepared to consider 
Martin’s comments. One of his objections is that, 
when the authors’ numerical procedure is applied to 
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that approach the solid line in Fig. 2 as the incre- 
ments in time and distance approach zero. 

The Buckley-Leverett flow considered by the 
authors is a very simple flow compared to what the 
partial differences can represent. There may be 
many discontinuities present in the porous media. 
They can arise within the body, and they can begin 
at the boundaries even after the flow has been in 
progress for some time. Since the authors have 
failed to present a completely satisfactory demon- 
stration that their methed yields the correct results 
for a very simple case, how can one be sure the 
method will be applicable in a very involved case 
which may arise in the analysis of a natural reser- 
voir? 
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the case of complete pressure maintenance, it does 
not approximate what he calls the ‘“‘true analytical 
solution’’ of the differential equation. The solution 
to which he refers was given by Buckley and Leverett 
and contains multiple values of saturation at some 
locations in space. Obviously, this is not a physi- 
cally significant solution of the problem and so 
must be rejected. The generalized solution which 
does occur in nature is formed by replacing the region 
of multiple-valued saturations by a saturation dis- 
continuity so located that the resulting solution 
satisfies the law of conservation of mass. The 
authors selected the difference analogue of the 
differential equation with the objective of approxi- 
mating the solution having physical significance. 
The comparison of Fig. 2 gives a measure of the 
degree of achievement of this objective. Martin is 
critical of the accuracy of the approximation in the 
vicinity of the saturation discontinuity. He states, 
‘*.,. It is not clear that the finite-difference method 
would yield results that approach the solid line in 
Fig. 2 as the increments in time and distance ap- 
proach zero”, The authors take issue with this state- 
ment. For a perfect solution, the length of theregion 
of the reservoir containing saturations between zero 
and the frontal saturation should be zero. The data 
of Fig. 2 show that when a distance increment of 
1,000 ft was used in the numerical procedure, this 
region occupied 22 per cent of the reservoir (from 
7,000 to 9,200 ft). For a distance increment of 500 
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ft, it occupied 14 per cent; for a distance increment 
of 250 ft, it occupied 8 per cent. Thus, halving the 
size of the distance increment while holding At/Ax 
constant approximately doubles the slope of the 
saturation curve in the vicinity of the front. This 
experimental relationship of error to grid size is the 
same as predicted by an elemental truncation error 
analysis of the difference equations used. This 
analysis shows that the size of the error should be 
proportional to Ax. The authors feel that this evi- 
dence makes clear the convergence of the numerical 
method as the time and distance increments approach 
zero. 

Martin’s final criticism is that the testing of the 
procedure was not sufficient to verify its accuracy 
in ‘‘a very involved case which may arise in the 
analysis of a natural reservoir’. While the testing 
of an experimental method is never complete, the 


validity of the present method has been demonstrated 
for a wide variety of conditions. To date, the authors 
have not encountered a case for which it fails. If 
evidence of such failure is obtained, the authors 
will welcome the information. 
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ABSTRACT 


The effect of a sand-filled vertical fracture of 
limited radial extent and finite capacity (fracture 
capacity is the product of the permeability and 
width of the fracture) on the flow behavior of a 
cylindrical reservoir producing an incompressible 
fluid through a centrally located well has been in- 
vestigated mathematically. The shape of the lines 
of equal pressure near the fracture is essentially 
independent of the size of the reservoir, provided 
that the field radius is of the order of the fracture 
length or larger. For reasonable values of produc- 
tion rates and of fluid, reservoir and fracture 
properties, the total pressure drop between the end 
of the fracture and the well is generally negligible 
compared with the pressure drop in the reservoir. 

With regard to production response, the effect of 
vertical fractures can be represented by the pro- 
duction response of an equivalent or effective well 
radius. For a high-capacity fracture, the effective 
well radius is a quarter of the total fracture length, 
decreasing with the fracture capacity. 

When invasion effects are simulated by decreas- 
ing the width of the damaged zone with distance 
from the well, the effect of formation damage around 
a fracture on the production response is not so 
serious as indicated by the literature. This sug- 
gests that frac fluids with a conventional filter-loss 
response are better than high-spurt-loss frac fluids, 
provided the effective permeability of the damaged 
zone is the same. 

INTRODUCTION 


This paper considers the effect of the fracture 
Capacity, as well as the formation damage which 
can result from fracture treatments, on the produc- 
tive capacity of vertically fractured wells. Other 
publications, notably those of van Poollen,!+2 con- 
sider these same effects. 

In addition to providing more general results for 
vertical fractures than are available from the liter- 
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ature, the present paper gives the equivalent well 
radius of fractures having different lengths and 
capacities and, also, includes pressure distribu- 
tions in and around the fractures. 

The effect of a damaged zone around a fracture 
on the production response was not found to be so 
great as that reported by van Poollen.? This differ- 
ence probably stems from the fact that we consider 
a damaged zone which is widest (but is still small) 
near the well and thins out toward the extremities 
of the fracture, whereas van Poollen considers a 
damaged zone having a uniform width for the entire 
fracture length. 

Simple, but adequate, equations which describe 
the effect of these variables on production response 
are presented (in Appendixes A and B). Thus, re- 
sults can easily be extended to values of the 
variables not specifically considered here. 


IDEALIZATION AND DESCRIPTION 
OF THE FRACTURED SYSTEM 


It is assumed that a horizontal oil-producing 
layer of constant thickness and of uniform porosity 
and permeability is bounded above and below by 
impermeable strata. The liquid is incompressible. 
It is assumed that, after the well is fractured, the 
cylindrical outer ‘‘boundary’’ is at a uniform po- 
tential, provided that it is not too near the fracture. 
The fracture system is represented by a single, 
plane, vertical fracture of limited radial extent, 
bounded by the impermeable matrix above and below 
the producing layer (reservoir).? 

Fig. 1 indicates the general three- dimensional 
geometry of the fractured reservoir described above. 
If gravity effects are neglected,.the flaw behavior 
in the reservoir will be independent of the vertical 
position in the oil sand. This means that the flow 
behavior in the fractured reservoir can be described 
by the two-dimensional flow behavior shown in the 
horizontal cross section of the reservoir, Fig. 2. 

PARAMETERS 

The effect of a fracture on the pressure distribu- 

tion in the fracture itself and in an otherwise un- 


damaged (no skin) reservoir indicates that the pres- 
sure distribution is a function of three parameters. 
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The first of these can be described as the ratio of 
the ability of the formation to carry fluids into the 
fracture to the ability of the fracture to carry fluids 
into the well, and it is defined by 

a= 

7 (formation permeability, k) (fracture length, 2L) , 





4 (fracture permeability, ky)(fracture width, w) 
or 
mk2L 
cee . 
4kyw 
The last two parameters describe only the geom- 
etry of the system, that is, the size of the wellbore 
and drainage radius relative to the extent of the 
fracture from the well axis (L). These parameters 
are denoted by rp and rep. * 


an tee «6 ee 


Thus, for a given geometry, the effectiveness of 
the fracture is described by the relative capacity 
parameter a. For a very effective fracture — one 
having a great ability to carry fluids to the well — 
the relative capacity parameter is very small and 
approaches the limiting value of zero for infinite 
fracture permeability. Likewise, for a very ineffec- 
tive fracture—one having a low capacity—the param- 
eter a is relatively large, approaching infinity for 
the limiting case of no fracture. 

Practical values of the parameter a appear to be 
less than 0.1 when the fracture permeability is the 
same as that of the sand used in the treatment. It 
is known, however, that fine particles from the res- 
ervoir, and in some cases other suspended solids 
in the fracturing fluid, will mix with the fracturing 
sand and will greatly reduce its permeability. In 
these cases, the values of the capacity parameter a 
could be much larger than 0.1. 

When the formation is damaged as a result of 
fracture treatment, the permeability damage and the 
depth of the damaged zone (relative to the fracture 
length) are added parameters. In this case, well 
effects and pressure drops in the fracture have been 
neglected to simplify the presentation of results. 





*All dimensionless distances are in terms of L, half the 
fracture length. For example, rep = re/L. 
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RESULTS AND DISCUSSION 


For formations with no skin, this study provides 
two types of information — one, such as the pressure 
distribution in and near the fracture, gives a general 
picture and understanding of the action of the frac- 
ture on the pressures behavior in the formation, and 
the other describes quantitatively the action of the 
fracture on such items as the magnitude of the pres- 
sure drop in the fracture and improvement in pro- 
duction rate. 

For formations damaged as a result of the frac- 
ture, production rates are compared to those which 
would have been obtained without damage. These 
results will be treated separately at the end of this 
section. 


PRESSURE DISTRIBUTION IN 
AND NEAR THE FRACTURE 

The dimensionless pressure drop between the 
end of the fracture and any point in the fracture is 
essentially independent of the size of the reservoir, 
i.e., of the length parameter rep. However, this 
pressure drop is strongly dependent on the value of 
the relative capacity parameter a. 

Fig. 3 shows the variation of this dimensionless 
pressure drop in the fracture with the dimensionless 
distance from the well, This dimensionless pres- 
sure drop is given by 
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TURE. 
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In Eq. 2, py — px is the pressure drop between 
the end of the fracture at x = L, in the co-ordinate 
system of Fig. 2, and any point in the fracture a 
distance x from the well axis. 

The curves presented in Fig. 3 indicate that as 
the fracture becomes less important (increasing a), 
the dimensionless-pressure-drop distribution in the 
fracture approaches that corresponding to strictly 
radial flow into a well of zero radius.* On the other 
hand, as the ability of the fracture to carry fluids 
increases (decreasing a), the pressure variation in 
the fracture becomes smaller, as indicated by the 
small values of AP; x. Thus, the pressure drop in 
the fracture is negligible unless the fracturing sand 
is contaminated with fine particles from the forma- 
tion, either during the course of the fracture treat- 
ment or after the well is put on production. 

It is interesting to note from Fig. 3 that, in general, 
the pressure gradients in the fracture at the well 
are finite even though the well has a zero radius. 
The only exception is for the limiting case of no 
fracture given by the curve for infinite a, for which 
case the pressure behaves as In r near the well. 
Thus, for a given production rate, the effect of the 
fractures is to reduce the pressure gradients near a 
well, 

Pressure distributions in the reservoir were de- 
termined for a few values of the capacity parameter 
aand for a value of the length parameter rep = 2.40. 
Values given at the outer boundary are numerically 
equal to AP. 4, which is discussed in the next sec- 
tion. These pressure distributions in the reservoir 
were obtained in a few hours by using relaxation 
methods,® except for values of a equal to zero and 
infinity, for which cases the shapes of the equi- 
pressure lines were known. 

Figs. 5, 6 and 7 show the pressure distribution 
in the reservoir for a=0.1, 1.0 and 10, respectively. 
The pressure distribution for the case of no pres- 
sure drop in the fracture (a = 0) is given by con- 
focal ellipses, as shown in Fig. 4. For infinite a, 
the fracture is of no influence, and the equipressure 
lines are circles concentric with the well axis (Fig. 
8). To give the pressure distribution for infinite a 
in a manner which could be compared with the curves 
in Figs. 4 through 7, that is, for rep = 2.40, the 
the equipressure lines in Fig. 8 were determined 
for a fictitious fracture length of 416 ft in a reser- 
voir of 500-ft radius drained by a well of 1/4-ft 
radius. In Figs. 4 through 7, the value of the equi- 
pressure lines represents the dimensionless pres- 
sure drop between the equipressure lines and the 





*The results of Fig. 3 indicate that the pressure distribution 
for a > 100 is very nearly the same as that for radial flow. When 
ever the permeability of the fracture is the same as that of the 
reservoir, the fracture is trivial because it is not effective. 
This means that the pressure distribution for such an ineffectual 
fracture must be that associated with values of a > 100. How- 
ever, for such fractures, the parameter a is given by a = 7L/2w, 
so that the method used in this study to obtain the pressure 
distribution appears to be limited to cases where the length/ 
width ratio of the fracture is larger than 130, This is not com 
sidered to be a practical limitation. 
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well of zero radius. In Fig. 8, the value of the equi- 
pressure lines represents the dimensionless pres- 
sure drop between the equipressure lines and the 
well of 1/4-ft radius. 
SPECIFIC PRODUCTIVITY INDEX 

It has been found that the pressure drop between 
the outer boundary and a well of radius r,, can be 
be represented as the difference of the total pres- 
sure drop to a well of zero radius, with a correction 
term added to represent the pressure drop between 
the well of radius 7, and one of zero radius. 


AP ow (4,7%ep. Twp) ™ APS 0 (4,7%ep) io 


AP, o (4%uyp)s - - « 
Here, APs og: the dimensionless pressure drop 
between the outer boundary and a well of zero radius, 
is a function of a and rep only; and AP, ,, the 
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dimensionless pressure drop between a well of 
radius r,, and a well of zero radius, is a function of 
a and ryp only. We note that the main contribution 
to the two terms in the right-hand side of Eq. 3 is 
in APe,o, which is plotted in Fig. 9. It should be 
noted that AP, , is also the reciprocal of the spe- 
cific productivity index SPI defined by 


a | 
kW be — Po) 


The term AP. plotted in Fig. 10 is only a cor- 
rection term which takes into account the effect of 
the well radius. 

For values of rp larger than two, the curves of 
AP e,o given in Fig. 9 can be represented very well 
by 1 
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where ry,p(a) =7,,/L is only a function of a. Values 
of the effective well radius rj,p(a) are given in 
Fig. 11, showing that for fractures of infinite ca- 
pacity, the effective well radius is always equal 
to a quarter of the fracture length, (2L)/4. Smaller 
values of the effective well radius correspond to 
less effective fractures. 

The well pressure used in the correction term 
APw,o given in Eq. 3 is approximated by the arith- 
metic average of the pressure along the fracture 
and the pressure along the normal to the fracture, 
with both these pressures being evaluated at a 
distance r from the center of the fracture. Fig. 10 
shows the correction term AP, , obtained as de 
scribed. We note that the correction term is unim- 
portant for reasonably effective fractures (a < 1), 
Figs. 6, 7 and 8 indicate that averaging the well 
pressure in the manner described is quite adequate 
for a2 1. 

MAGNITUDE OF THE PRESSURE 
DROP IN THE FRACTURE 

It is common practice now to consider that the 
pressure drop in the fracture is small compared to 
the pressure drop in the field. In fact, it is generally 
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assumed that there is no pressure drop in the frac- 
ture.4°5 It is desirable to determine the conditions 
of fracture effectiveness and length that will permit 
neglecting the pressure drop in the fracture (say, 
less than 10 per cent of the pressure drop in the 
field). A plot of the ratio of pressure drop in the 
fracture to that in the field is given in Fig. 12 for 
for a well of zero radius. Consideration of the small 
correction terms for the well radius would only com- 
plicate the presentation of the results without 
changing their character. The ratio of the pressure 
drops plotted in Fig. 12 as a function of a and rep 
has been obtained from the results presented in 
Figs. 3 and 9. The results presented in Fig. 12 
show that the pressure drop in the fracture is 
negligible (less than 10 per cent of the pressure 
drop in the field) for values of a <_0.1 and for 
fracture lengths smaller than the field radius. 
Since most of today’s fractures probably satisfy 
these conditions, it appears that pressure drops in 
clean fractures can be considered negligible in most 
present cases. 

However, the curves in Fig. 12 show that the 
ratio of the pressure drops is strongly dependent on 
the values of a and rep. Thus, for values of a = 0.5 
and fractures as short as 1/40 of the field radius, 
the pressure drop in the fracture is between 10 and 
30 per cent of the pressure drop in the field. This 
means, of course, that the pressure drop in fractures 
can become an appreciable fraction of the pressure 
drop in the field if the propping agent becomes con- 
taminated with fine particles. Also, the present 
trend in using larger volumes of fracture treatment 
(as well as increasing the sand content) will probably 
result in larger values of a and smaller values of 
tep. It follows that appreciable pressure drops in 
the fracture may become more frequent in the future, 
even when the propping agent remains clean. 


IMPROVEMENT IN PRODUCTION RATE 
The improvement in production rate q/q is given 
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EFFECTIVE WELL RADIUS ryp 


As an example, the improvement in production rate 

is determined below for a well of 1/4-ft radius in 

the middle of a 200-ft fracture of infinite capacity 

(a = 0) and ‘‘draining’’ from a radius of 500 ft. Us- 

ing the results of Figs.9 and 10 in Eq. 5, we obtain 
1 in 2000 


7 


Y . = 33. 

q 0.365 (from Fig. 9) — 0(from Fig. 10) 
Whenever the effect of the well radius is small, the 
improvement in production rate can be obtained 
most simply by using the effective well radius of 
the fracture in the formula 
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For the conditions described above (a = 0), 
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as was also obtained by the more general method 
given by Eq. 5. 

This means that the maximum increase in pro- 
duction rate per unit pressure drop for such condi- 
tions cannot be expected to be more than 3.3 times 
the rate obtained before fracturing, assuming that 
the well is not damaged before the fracture treat- 
ment. a 

Increases in the production rate (SPI) can gener- 
ally be associated with increases in the effective 
radius of the well, as can also be obtained, say, by 
acid treatments. Whenever the formation is homog- 
eneous and is not fractured during the acid treat- 
ment, the acid enters the formation radially, enlarging 
the well radius equally in all directions to a new 
effective well radius of value ry. The example 
given above suggests that, if a well is stimulated 
by acidizing the formation, the acid treatment must 
result in an effective radius of 50 ft before the 
improvement of production rate per unit pressure 
drop can match the improvement due to a single 
fracture of infinite capacity having a total length of 
200 ft. 

A similar interpretation can be placed on the 
other values of effective well radius associated 
with the values of the parameter a given in Fig. 9. 
These values of ry,p indicate that, for large values 
of a, acidizing shows more promise than fracturing 
and that, for relatively small values of a, long frac- 
tures show more promise than acid treatment. For 
intermediate values of a, the relative merits of the 
two methods, according to the interpretation given, 
will largely depend on the size of the fracture and 
of the acid treatment. 

Still another interpretation can be given to the 
effective well radius. It is interesting to note that, 
for a hypothetical case of an actual well radius of 
100 ft (corresponding to a large number of distinct 
finite vertical fractures, each 200 ft in length, cen- 
tered at the well axis), the maximum improvement 
in productivity index is 4.7 times the productivity 
index given by a well of 1/4-ft radius draining a 
field of 500-ft radius. Since the improvement due to 
a single fracture for the same conditions has been 
shown to be 3.3, this factor of 4.7 indicates that, 
after one long fracture is created, the additional 
gain by multiple radial fractures (assuming that it 
is possible to form such a system of fractures) will 
be small. This, in turn, suggests that one long frac- 
ture may be preferred to several shorter ones. 


MAXIMUM PRODUCTION RATE FOR 
A GIVEN FRACTURE VOLUME 
For a given fracture volume, an increase in the 
width results in a shorter fracture. These changes 
have opposite effects on the production rate, as 
can be determined from Fig. 9. This raises the 
question of which parameter values result in a max- 
imum rate of production for a given fracture volume. 
This question can be answered by finding the 
value of a at which q at constant fracture volume 
(constant ar2p) has a maximum. The result is that, 
for a fracture of any volume and a production well 
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of zero radius, the maximum production rate is ob- 
tained when a is about 1.25. Thus, from the defini- 
tion of a, the optimum length- width ratio of any 
fracture is 


‘ 
a Jee 
w k 


Thus, for a fracture volume of 6.4 cu ft/ft of net 
pay, when ky = 100 darcies and k = 0.01 darcy, 
the optimum dimensions for this fracture volume 
are w = 0.02 and 2L = 320 ft. Less reasonable 
values of w and 2L could be obtained with different 
(and yet practical) values of k and k;. 

If we could learn how to pant the average 
separation between fracture faces, fracture treat- 
ments could be planned to give optimum results 
not only on the basis of fracture area alone, but 
also upon consideration of the length-width ratio 
of the fracture. These results are valid only in the 
area below and to the left of the line labeled rep = 
2 in Fig. 13. Dashed curves have been extrapolated 
from rep = 2 to rep = 1,* but this portion of the 
curves cannot be used with confidence. 


EFFECT OF SKIN ON PRODUCTION RATE 

Results presented thus far consider no damage 
to the formation during the fracturing operations, 
This section will consider the possible loss of 
potential production due to a damaged zone around 
the fracture in a region of relatively low permea- 
bility, resulting, say, from the plugging of the 
formation by certain fluid-loss additives which may 
be used in fracturing treatments. 7 

To obtain an idea of the degree of impairment in 
production rate due to such a damaged zone, we 
consider only the case in which the fracture capac- 
ity is infinite (a = 0). Physically, the decrease in 
production rate with the damage depends on the 
flow resistance of the skin compared with that of 
the formation. It must be noted that the flow resist- 
ance depends on the flow pattern. In an unfractured 
well, this flow resistance is largest nearest the 
well because the fluids converge radially into the 
well over its relatively small surface area, Any 
damage to this zone of high resistance may thus 
result in drastic reductions in production rate. In 
a fractured well, the flow is distributed over the 
large surface area of the fracture, resulting in al- 
most linear flow in the formation near the fracture 
and, consequently, a small flow resistance. In the 
latter case, damage to the formation near the frac- 
ture may not change the total flow resistance appre- 
ciably and, thus, could not be expected to reduce 
the production rate significantly. Details of the 
solution of this problem are found in Appendix B. 
A horizontal cross section of the reservoir-fracture- 
skin system is shown in Fig. 14. 

The ratio of production rate after damage to un- 
damaged production rate is 





*For rep = |, the fracture extends to the (pre-fracture) draim 
age radius of the field, 
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as developed in Appendix B. 

Figs. 15 and 16 show the calculated ratio of 
production rate after damage to undamaged produc- 
tion rate for fracture lengths of 200 and 40 ft, 
respectively, in a reservoir of 500-ft radius. A 
permeability reduction of 99 per cent in a damaged 
zone of l-in. average depth results in a reduction 
of only 4 and 13 per cent for the 200- and 40-ft 
fractures, respectively. Within this range of fracture 
length, permeability reduction and depth of damaged 
zone, the effect of the skin on the production rate 
is relatively unimportant. However, for a similar 
damaged zone, a reduction in permeability of 99.9 
per cent results in a reduction in production rate of 
32 and 58 per cent for the 200- and 40-ft fractures, 
respectively. 

Fig. 15 also gives one of the data points of Ref. 
2, which shows more production impairment than 
our results show. This difference is believed to 
result from the assumed shape of the damaged zone. 
In this paper, the skin is widest at the well and 
thins out toward the ends of the fracture. This 
simulates the damaged zone due to filtrate invasion 
which could be obtained with frac fluids having 
conventional filter-loss response (that is, the loss 
is proportional to the square root of the time). In 
Ref. 2, the width of the skin is taken to be constant 
over the extent of the fracture. This simulates the 
damaged zone due to filtrate invasion which could 
be obtained with frac fluids having an initial spurt 
loss and a small subsequent loss. In this latter 
case, the depth of the damaged zone and the degree 
of damage would be approximately uniform over the 
surface area of the fracture. On the basis of pro- 
ductivity impairment, the use of the frac fluid hav- 
ing a conventional filter-loss response is indicated. 


Vg= Vy k/ Kg 


Vg = Ol Vpk/ Ky 


= 0,01 vyk/ky 


SPl, AT CONSTANT FRACTURE VOLUME 





0.1 1. 
RELATIVE CAPACITY PARAMETER, a 


FIG. 13— CURVES SHOWING THE OPTIMUM PRODUC- 
TION RATE. 





Results of a similar nature were also determined 
for unfractured wells surrounded by a damaged zone. 
Fig. 17 shows the calculated ratio of production 
rate after damage to undamaged production rate for 
an oil sand which has been in contact with fractur- 
ing fluids but which has not been fractured. 
Muskat’s 8 equation (7.3.11) was re-arranged as 


Vfd _ 1 
qy In(1+ pa 
Tw 


1 Oe ‘- i) ———_£- 
ka In “¢ 
Tw 


(9) 





to obtain these results. For the case of a well of 
3-in. radius having a damaged zone 2-in. deep, a 
permeability reduction of only 90 per cent reduces 
the production rate by 37 per cent. Thus, we can 
conclude that the production rate from any sand 
section which fails to fracture will be reduced 
severely, especially if the sand section is separated 
from the fractured sections by shale breaks. This 
loss in potential production may be obscured by 
the increased productivity of the fractured sections. 


CONCLUSIONS 


Some general conclusions drawn or suggested by 
the results of this study are briefly discussed here. 
It must be kept in mind that these conclusions are 
valid only within the limitations of the idealizations 
used in this investigation. However, they serve at 
least as valuable guides in other situations. 

A high-capacity vertical fracture has an effective 
well radius of one-fourth the fracture length. This 
equivalence can be used to advantage to show that 
the production rate resulting from a 200-ft, high- 
capacity fracture in a field of 500-ft radius drained 
by a well of 1/4-ft radius is 3.3 times the produc- 
tion rate before fracturing. Changing either the well 
or field radius or the fracture length by a factor of 
two results in production rates from 2.5 to 4.7 times 
the rate before fracturing. The order of magnitude 
of this improvement in production rate is believed 
to correspond to that obtained at semi-steady state 
in depletion-type reservoirs. The equivalent well 
radius for fractures of lower capacity is given in 
the paper. 

The improvement in production rate due to a high- 
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FIG. 16 —-LOSS OF PRODUCTION DUE TO DAMAGE AROUND A SHORT VERTICAL FRACTURE. 
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capacity fracture is equal to the improvement in 
production rate due to an acid treatment which re- 
sults in an effective well radius of one-fourth the 
fracture length. This suggests that fracturing is to 
be preferred to acidizing if the formation of long 
fractures appears probable. This conclusion is 
based on the premise that acid is injected radially 
into a homogeneous formation (no natural fractures 
exist) and that the formation is not fractured during 
the acid treatment. 

The pressure drop in the fracture is small (less 
than 10 per cent of the pressure drop in the field) in 
most present cases, that is, for a< 0.1 and fracture 
lengths less than the field radius. However, pres- 
sure drops in the fracture of the order of 25 percent 
(and more) of the pressure drop in the field may be 
expected for values of a = 0.5. Such values of a 
may be realized if fine particles mix with the prop- 
ping agent or possibly if the present trend in using 
larger volumes of fracture treatment and sand con- 
tent continues. 

Extensive damage to the formation around a frac- 
ture (99 per cent permeability reduction for a depth 
of 1 in.) may not necessarily result in appreciably 
poorer production rate response (up to 13 per cent 
for rp /L = 5). When the depth of the damaged zone 
is proportional to the filtrate invasion, the de- 
crease in production rate resulting from a frac 
fluid with a conventional filter-loss response ap- 
pears to be less than that resulting from a high-spurt- 


loss frac fluid. However, the production rate from any 
producing section which fails to fracture may be 
reduced appreciably (25 per cent for a 90-per-cent 
permeability reduction for a depth of 1 in. around a 
3-in. well radius), Thus, in fracturing a long pro- 
ducing section, loss in potential production due to 
damaged zones in sections which fail to fracture 
may be obscured by the increased productivity of 
the fractured sections. 


NOMENCLATURE 
a = (n/4) (2Lk/wk ), dimensionless 
hb = net formation thickness, cm 
k =effective formation permeability to oil, 
darcy 
L =half the fracture length, ft 


In = logarithm to the base e 


= pressure, atm 


p 

q = production rate, cc/sec 

r = distance from well axis, ft 
, 
w 


= effective well radius of a fracture, ft 


~” 
'U 
— 
" 


qu/khAp, dimensionless specific produc- 
tivity index 

V; = volume of the fracture, cu ft 

Vp», = bulk volume of the formation, cu ft 
w = average width of the fracture, ft 

Wq = average depth of the damaged zone, ft 


x = distance from well, ft 
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AP;,; = kbp; — 2; )/qu, dimensionless 
Ap = pressure drop between outer boundary and 
well, atm 
pb; — pj =pressure drop between two points sym- 
bolized by i and j, atm 

# = viscosity of produced fluid, cp 

SUBSCRIPTS 

D =dimensionless, all distances are referred 
to L, one-half the fracture length 

d = damaged 

e = external boundary 

f = fracture 

L =end of fracture 

O = well axis 

w = wellbore 

x = distance along fracture 
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APPENDIX A 


METHOD OF SOLUTION USED TO OBTAIN 
THE PRESSURE DISTRIBUTION IN A 
VERTICALLY FRACTURED SYSTEM 


As a consequence of the assumptions made in 
the text with regard to the fluids and the formation, 
the pressure obeys Laplace’s equation both in the 
fracture and in the reservoir. The detailed mathe- 
matical procedures used in solving Laplace’s equa- 
tion to obtain the pressure distribution in the for- 
mation and in the fracture are given herein, but it 
seems desirable to indicate first two pertinent 
points about the method of solution. 

One of these points is that the circular outer 
boundary is replaced by an ellipse having an almost 
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circular shape. The ratio of axes of the ellipses 
chosen depends on the fracture length but is in all 
cases greater than 0.88. It was actually shown 
that, for rep 2 2, the use of the slightly elliptical 
boundary has a negligible effect on the value and 
shape of the interior equipressure lines given in 
Figs. 4 through 7. 

The second point is that the reservoir is divided 
into two regions of different permeability, Regions 
I and II, as shown in Fig. 18. Region I represents 
the fracture of permeability ky, and Region II repre- 
sents the formation of permeability k. 

Initially, Region I is represented by the interior 
of an ellipse with foci at x = +L. The circular 
outer boundary of radius re is replaced by an almost 
circular ellipse of equivalent radius, also having 
foci at x = +L. Region II, between the outer bound- 
ary and the inner ellipse, represents the formation 
proper. The pressure distribution in these two 
regions is then obtained by making use of conformal- 
mapping transformations and Fourier series. In the 
final stage of the method of solution, the minor 
axis of the inner ellipse is reduced to zero in such 
a way that the average width of the ellipse times 
its permeability corresponds to the fracture capac- 
ity. The inner ellipse of zero minor axis degenerates 
to the straight line with ends at x = +L and, thus, 
represents the fracture shown in Fig. 2. This is 
the method used to obtain the pressure distribution 
for the fractured system shown in Fig. 2. A some- 
what more detailed description follows. 


The pressure in Regions I and II obeys Laplace’s 
equation: 

7p, =O in RegioI ,........ (Arl) 

V7p. =O in Region], ........ (A=2) 


with refraction conditions 
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FIG. 18—AUXILIARY MODEL OF RESERVOIR SYSTEM. 
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‘ 1 = op2 - (A-3) 

! on an 
on the boundary of inner ellipse, where derivatives 
are taken normal to this boundary, and boundary 


conditions 


1 _ 9 


on y=0,x%40..... (A-4) 


oo. 2. 988s, s.c40-ss0 Ce 


2|F =| 


Pz =C, (a constant) on the outer boundary, and q/b 
is the production rate per unit thickness at the well 
located at x = 0, y =0. 

It is difficult to apply these conditions in a 
Cartesian co-ordinate system on elliptical curves; 
thus, the transformation z = x + jy = L cosh €#=L 
cosh (& + jn) is used to map the two regions in the 
first quadrant of the reservoir shown in Fig. 18 into 
the rectangles shown in Fig. 19. 

The pressure in the rectangular Regions | and II 
of Fig. 19 still obeys Laplace’s equation and the 
refraction conditions at the common boundary be- 
tween the two regions, represented now by & = &,. 
The boundary conditions become 


Op1 
an 


and 


=0 on n=0, 1= > . . (A-6) 


3|$ 


+ 


_* on €=0 ..... (A*7) 


because of symmetry, and 


| 


ut 
Py =Cg on Ex€., where,;~ = 2e ©. (A-8) 


q/b is the production rate per unit thickness at the 
well located at €=0, n = 2/2. 

To decrease the likelihood of obtaining a slowly 
convergent series for the pressure distribution ow- 
ing to the singular behavior of the production well 
(point sink) at D (€ =0, n = 7/2) in Fig. 19, it was 
considered desirable to define reduced pressures 
by subtracting a solution having the same singular 
behavior at D and satisfying the boundary condi- 
tions but not the refraction conditions. Thus, the 
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reduced pressures are 
b; (€,n) = kyo (Em) - pg (Em)... « (A-9) 
2 (E,n) = kpa(E,n) -— py (En) . . . .(A-10) 


where p, (&,n), the pressure distribution having the 
desired singular behavior at D and satisfying the 
boundary conditions, is given by 


K 
Ps (€,m) = —q, Re ln dn (Sf 
e 





m) , .(A-11) 


where q, is the production rate in reduced units 
(9, = qp/2mb), and dn is one of the Jacobian elliptic 
functions of parameter m. The parameter m is chosen 
such that 


7 K’(m) 


4, K(m)’ 
where Km) and K(m), related by K“(1 — m) = K(m), 





- (A-12) 


are the complete elliptic integrals of the first kind. 
Properties of these functions can be found in any 
standard text on the subject, and the method for 
obtaining this solution is given in Ref. 9. 
The boundary condition at € = &, requires that 
Co = 0s (En) = -te In (1 -m) ; 
but, although p, (€,7) satisfies all the boundary con- 
ditions of p, and pp, it does not satisfy the refraction 
conditions. Therefore, the boundary and refraction 
conditions for the reduced pressures are 


«no n=0,9= + (A-13) 
on on : 
Ody 
a a =0.... (A-14) 
z on g ( 
gy =0 on f=€- + + + (A-15) 


kd, = kyo + (ky —k) Ps 


OS. 28 
aé dé 


on &= 6. (A-16) 


The use of Fourier series is particularly well 
suited to obtain the reduced pressure distributions 
, and do. These are given by 


$1 (&y) = _ Ay cos 2nn cosh 2n€ , . (A-17) 


$2 (&n) = > B,, cos 2nn sinh 2n(&--4 
opty: hike ae yt (A-18) 
where these expressions satisfy all conditions 
except the refraction conditions. By making use of 
the refraction conditions given by Eq. A-16, it is 
possible to determine the constants A, and B,, so 
that 
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$1 (Sm) = (ky - —k) lead In(1 —m) + 


kS 
n=1 [k +k, tanh 2n(€,- 1) tanh 2n€,] cosh 2néy 





K,,(0) cos 2nn cosh 2n€ 


- (A-19) 


and 
$y (€,n) = (ky a k)qo - 
2 K,,(0) cos 2nn sinh 2n(E-&) tanh 2né; 
ant [k+k, tanh 2n( &-€,) tanh 2n€, J cosh 2 Ee ) 








. (A-20) 
with 
gig TE a. 
n n ml 1+(—4,)" 1- q? 1? 
. (A-21) 


and where q; = e~*¢e = (L/2r, 4 is very small. 
From Eqs. A-9 and A-10, the pressure in the 
fracture (p;) and in the formation (p ) are given by 


ai _ bs (Sn) a ay. 
p, (én) ie (p, + $4) <a (a “ i) 


Yd Rasnt-oaie 
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ee K,,(0) cos 2nn cosh 2n€ | 
=1 [k+k, tanh 2n ,—€,) tanh 2n € | cosh 2n f 





- (A-22) 
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p, (én) 
p(n) => (b, + $2) = 





x ©K,(0) cos 2nn sinh 2n (€,-€) tanh 2né, 
rc [k+k, tanh 27K €-&) tanh 2n& ]cosh 2n &-€,) 


. (A-23) 


These solutions given for the pressures in the 
formation and in the fracture are valid for any frac- 
ture of elliptical shape as shown in Fig. 18. In 
practice, however, the fracture is narrow and long 
and has a high permeability. For these conditions, 
the reservoir pressure p(&) is given by 


p (En) = Limit po (Em &2,6,k, ky) » (A-24) 
€; +0 7 
ky > co 
2L Ek, = Ga’ 











_bs(Sn) , 
k 
24. 9 _ncos 2nn sinh 2n (& - €) K.(0) 
k n=1 (a+ 2ntanh 2n&) cosh 2n€, ” °’ 
a - (A-25) 
where 
2Lk 
- 2k .(A=26) 


and Ca’ is the permeability of the fracture times its 
maximum width, the minor axis of the ellipse. 

In a similar manner, the pressure in the narrow 
fracture can be found to be 


bm) = “fe [f+ += In (1 —m) + 


cos 2nn 
1 a + 2n tanh 2n€, 





oo 
a> Kq(0)] (A-27) 
n= 
Although the fracture capacity Ca’has been defined 
above in terms of the maximum width of the fracture, 
it may be more reasonable to define a new fracture 
Capacity Ca in terms of the average fracture width 
w. Since the average width of such a fracture is 
7/4 times its maximum width, it follows that the 
average capacity of the fracture Ca is given by 


Ca = Ca’ SP narig ne! Jeune de aenay ce te Renee 
and 
— z 2Lk .(A-29) 
wks 


Eqs. A-25, A-27 and A-29 were used for evaluat- 
ing the results discussed in the text. 


APPENDIX B 
EFFECT OF A SKIN ON PRODUCTION RATE 


Use is made of the methods discussed in Appendix 
A to investigate the effect of a skin on the produc- 
tion rate. For the case of infinite fracture capacity 
(a =0), the entire fracture is at a uniform pressure 
equal to the well pressure. We shall call this frac- 
ture pressure fy. The damaged zone around the 
fracture shown in Fig. 14 and the rest of the reser- 
voir can be transformed into rectangular sections 
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as shown in Fig. 20. 

We note from Fig. 20 that the pressure distribution 
in & should be the same for all values of 7; that is, 
the pressure distribution is only a function of é. 
This is because the pressure obeys Laplace’s equa- 
tion in both the &,n plane and in the x, y plane; 
therefore, the pressure distribution in the skin of 


permeability ky is 


pPa=zbr+ec€, 


and the pressure distribution in the reservoir is 


p=d+eé€ 


The constants b, c, d and e appearing in these 
equations are to be determined from the boundary 
conditions 


Pa = py at & = 0 (the fracture), 
p =p, at €=&,(the outer boundary of the field), 


and the refraction conditions 


Pa=P 
d, ; 

ka “Pd =k dp on &= &> (the boundary between 
dé dé the skin and the forma- 


tion). 


The quantity &; is related to the average depth 
of the skin by 





where 2L is the total length of the fracture. The 
quantity &, is related to the average field radius r 
by 


e 


The application of the boundary, refraction and Darcy 
equations leads to a production rate after damage 
given by 








_ 2nkb (Pe — by) 
y. 
oe 


p, 

( k 1) 4wq 
L ka aL 
Before the formation is damaged, ky = k, and the 
production rate is 


In 


2mkh (pe — py) 
L 
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Thus, the ratio of damaged to undamaged production 
rate is given by 








Vd _ I 
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DISCUSSION 


F. E. CAMPION 
MEMBER AIME 


The results of Prats’ analysis of flow behavior 
in vertical fractures are enlightening, particularly 
in regard to the effects of formation damage. The 
conclusions he has reached concerning the effects 
of damage to fractured wells appear valid within 
the framework of the assumptions set forth; however, 
further comments seem appropriate regarding the 
projection of these results to a physical prototype. 

The author has demonstrated resourcefulness in 
the treatment and handling of this problem, but we 
feel that the conditions assumed for evaluation of 
fracture damage are optimistic. There are a number 
of conditions that could render his results less 
favorable. For this reason, a direct extrapolation 
of his conclusions to a practical field situation 
may be difficult if the situation does not comply 
with the assumptions prescribed for the mathematical 
model, Frac fluids that damage formation permea- 
bility may not necessarily reduce production-rate 
response, but it is not difficult to envision a situa- 
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tion where they could. 

It appears from this study that the shape of the 
damaged zone surrounding the fracture is a critical 
factor in determining the fluid-flow behavior in the 
system and that much depends on the fluid-loss 
characteristics of the fracturing fluid. Since high 
spurt loss is not uncommon for frac fluids, less 
favorable production-rate response than that cited 
in the paper may occur as a result of the more uni- 
form distribution of damage along the fracture length. 

Considerations other than fluid-loss character- 
istics are also important in determining the suita- 
bility of fracturing fluids. The depth of formation 
damage may also be significant. Although depth of 
damage is limited, in many instances to 1 or 2 in., 
this is not consistent for all frac fluids. We have 
found wide variations in the depth as well as degree 
of damage and feel that assignment of representa- 
tive values for either of these factors is restrictive. 

The assumption of infinite fracture capacity 
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employed in the study may also limit the scope of 
the results. van Poollen! has shown that the effects 
of formation damage are compounded by reduced 
fracture capacity. Our experience shows that frac 
fluids which damage formation permeabifity may 
also reduce permeability of propping agents in the 
fracture. Although it has become common practice 
te neglect pressure drop in the fracture, this as- 





lvan Poollen, H. K.: ‘Productivity versus Permeability 
Damage in Hydraulically Produced Fractures’’, Drill. and Prod. 
Prac., API (1957) 103. 





Although the conditions assumed in the report 
for evaluating fracture damage may be optimistic 
in some cases, it is ciear that the conditions 
assumed in the other published works on this 
subject,2 namely that of a uniformly thick damaged 
zone, are definitely pessimistic. Even with the 
assumption of a uniformly thick damage zone 
over the entire fracture surface, results in Ref. 2 
indicate that the main contribution to the production 
is derived from the portion of the fracture surface 
which is farther from the well. When one considers 
that the time of contact between the frac fluid and 
the fracture surface is longest nearest the well, it 
is almost certain that the extent of damage and the 
depth of the damaged zone would be smallest far 
from the well, regardless of the type and properties 
of the frac fluid used. Thus, it is entirely reasonable 
to anticipate that the contribution to production from 
the sections of the fracture surface farther from the 
well under conditions existing in the formation after 
a frac treatment would be larger than indicated in 
Ref. 2. This means that the results of Ref. 2 regard- 
ing the effect of a damage zone on loss of production 
are pessimistic. 





AUTHORS’ REPLY TO S. E. CAMPION 


sumption is not wholly justified and conclusions 
based upon it tend to be optimistic. 

Laboratory measurements of permeability damage 
by fracturing fluids can do much to aid in the selec- 
tion of suitable frac fluids. Although these meas- 
urements may not be useful in a strict quantitative 
sense, they can serve as a basis for comparison of 
fracture fluids. Through such tests, the use of 
fluids which do extensive damage can be avoided 
and the risk of losing significant quantities of oil 
over a period of years can be minimized. 


It should be stressed that, although the general 
shape of the damaged zone can be surmised to thin 
out with increasing distance from the well, its 
exact shape in situ will probably never be determined. 
In fact, the shape of the damaged zone and the 
extent of damage will depend not only on the prop- 
erties of the frac fluid, but also on those of the 
reservoir rock fluids as well as on the pressure, 
saturation and temperature conditions existing 
in the reservoir at the time of the fracture treatment. 
For example, a gas saturation increasing towards 
the well could result in a damaged zone which 
thins out away from the well. Thus, it cannot be 
said that the results presented in the paper are 
generally optimistic, although they might be in 
some cases. By contrast, the results of Ref. 2 are 
almost certainly pessimistic. 


Although the paper does not consider the multiple 
effects of a damage zone and reduced fracture 
capacity, similar arguments suggest that the results 
presented for such conditions in the only other 
reference to the subject are also somewhat pessi- 
mistic. xe 
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